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Dynamics  Division,  Air  Force  Flight  Dynamics  Laboratory, 
Wright-Patterson  Air  Force  Base,  Ohio  under  Contract 
F33615-69-C-1339 .  This  research  is  part  of  a  continuing 
dynamic  and  aeroelastic  phenomena  for  rotor/propeller  powered 
V/STOL  flight  vehicles  under  the  Air  Force  Systems  Command's 
exploratory  development  program.  The  Project  Number  is  1370, 
"Dynamic  Problems  in  Military  Flight  Vehicles",  Task  Number  is 
137005,  "Prediction  and  Control  of  Flight  Vehicle  Vibration". 
Mr.  A.  R.  Basso  of  the  Aerospace  Dynamics  Branch  was  the 
Project  Engineer. 
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predicting  folding  proprotor  (FPR)  aircraft  dynamic  behavior, 
correlation  of  these  methods  with  model  test  data,  and  the 
results  of  the  dynamic  analysis  of  a  representative  FPR  design. 
The  second  volume  is  a  guide  to  the  digital  computer  programs 
and  contains  input  and  output  formats  and  FORTRAN  listings  of 
the  programs. 

Mr.  Troy  Gaffey  was  the  Bell  Helicopter  Company  Project 
Engineer.  The  authors  acknowledge  the  assistance  of  Dr.  Jing 
Yen  and  Mr.  Gottfried  Weber  in  providing  the  equations  of 
motion  and  in  checking  out  the  programs. 
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1971  for  publication  as  an  AFFDL  Technical  Report. 
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ABSTRACT 


This  report  describes  the  results  of  a  study  of  the  dynamic 
stability  and  response  during  the  feathering  and  folding 
of  the  blades  of  a  folding-proprotor  VTOL.  This  study  involved 
the  development  of  analytical  methods  of  predicting  the 
dynamic  characteristics  during  feathering  and  folding,  corre¬ 
lation  of  the  theoretical  methods  with  experimental  data,  a 
dynamic  analysis  of  a  representative  folding-proprotor  VTOL 
design,  and  a  parametric  study  to  identify  design  factors 
which  can  be  used  to  control  dynamic  characteristics  during 
feathering  and  folding.  Correlation  of  the  theory  with 
measured  dynamic  stability  and  response  characteristics 
during  feathering  and  folding  is  good,  and  indicates  the 
analytical  methods  can  be  used  with  confidence.  The  results 
of  the  dynamic  analysis  indicate  satisfactory  stability  and 
response  characteristics  can  be  achieved  in  a  66,000-pound 
gross-weight  class,  folding-proprotor  VTOL.  The  results  of 
the  parametric  study  are  summari2ed  in  terms  of  design  guide¬ 
lines.  This  volume  is  a  guide  to  the  computer  programs,  and 
contains  FORTRAN  listings.  Volume  I  contains  the  develop¬ 
ment  of  the  theory,  correlation  of  theory  with  experimental 
data,  and  the  dynamic  response  and  parametric  study. 
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SECTION  I 


INTRODUCTION 


This  volume  provides  a  guide  to  the  use  ot  the  computer  pro¬ 
grams  developed  during  this  sLudy.  It  also  contains  the  infor- 
mation  necessary  to  inodity  the  programs  or  to  adapt  them  to 
a  given  computer.  A  FORTRAN  listing  ot  the  programs  is  in¬ 
cluded  in  this  volume.  It  should  be  noted  that  the  program 
listings  arc  the  only  source  tor  the  complete  equations  of 
motion.  Section  II  of  Volume  I  provides  an  overview  of  the 
theory  mechanized  in  the  programs. 

Three  computer  programs  were  developed:  ARAP06 ,  which  treats 
that  portion  of  the  FPR  transition  mode  where  the  blades  are 
rotating,  and  DFAL17  and  DFAL18,  which  treat  that  portion 
where  the  proprotor  has  been  stopped  and  the  blades  are  being 
folded.  ARAP06  is  all  inclusive  in  that  it  may  be  used  to 
calculate  the  stability,  loads,  and  response  in  proprotor  mode 
and  during  blade  feathering.  DFAL17  is  used  to  calculate  the 
system  stability  during  folding;  DFAL18,  the  loads  and  response 
during  folding. 

Sample  cases  for  each  program  are  contained  in  Appendix  I. 

These  include  card  data  and  condensed  output. 
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SECTION  II 


USER  GUIDE  TO  PROGRAMS 


In  preparing  this  section  it  is  assumed  that  the  reader  has 
read  Section  II  of  Volume  I  and  is  familiar  with  the  math 
model  and  capabilities  and  limitations  of  each  program. 

This  user's  guide  consists  of  four  parts  for  each  program: 

(1)  an  input  format  identifying  the  information  contained  on 
each  card,  (2)  an  input  format  guide  explaining  the  function 
and  source  of  each  input,  (3)  an  output  format  and  guide, 
and  (4)  suggested  techniques  for  using  the  programs. 

A.  ARAP06  -  FPR  FEATHERING  DYNAMICS  ANALYSIS 

The  inputs  to  ARAP06  consist  of  two  types  of  data:  (I)  math 
model  configuration  and  control  options  and  (2)  vehicle 
descriptive  data.  The  options  allow  the  user  to  select  the 
math  model,  the  type  of  excitation  and/or  the  maneuver  to  be 
considered,  and  the  output  data.  The  descriptive  data  consist 
of  the  geometry,  certain  spring  rates,  mass  distributions,  and 
the  coupled  normal  modes  of  the  subject  vehicle. 

To  generate  the  normal  mode  input  data  the  user  must  have 
access  to  a  coupled  beam-chord-torsion  vibration  analysis  (to 
determine  the  normal  modes  of  the  wing/pylon  system)  and  to 
a  rotating  beam  vibration  analysis  (to  determine  the  blade 
normal  modes). 

Because  of  the  various  options  with  regard  to  the  math  model,  all 
of  the  input  data  will  generally  not  be  required.  For  example, 
in  many  cases  the  equations  representing  blade  elastic  motion 
may  be  deleted  without  significantly  influencing  the  results. 

In  those  cases  blade  normal  mode  data  are  not  required  as  a 
program  input.  Careful  selection  of  the  options  for  a  given 
case  can  greatly  reduce  the  data  preparation  effort. 

Figure  1  illustrates  the  card  order  in  a  ARAP06  data  deck. 
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Figure  1.  Card  Order  in  ARAP06  Data  Deck 
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1.  Input  Format 


CARD  OL  Title  Card;  Column  1-80,  alphameric. 


CARD  02 


CARD  03 


Format 

1615: 

Col. 

1 

— 

5 

IT1 

Col. 

6 

- 

10 

IT2 

Col. 

11 

- 

15 

IT3 

Col. 

16 

- 

20 

IT4 

Col. 

21 

- 

25 

IT5 

Col. 

26 

- 

30 

IT6 

Col. 

31 

- 

35 

IT7 

Col. 

36 

- 

40 

IT8 

Col. 

41 

- 

45 

IT9 

Col. 

46 

- 

50 

IT10 

Col. 

51 

- 

55 

IT11 

Col. 

56 

- 

60 

IT12 

Col. 

61 

- 

65 

IT13 

Col. 

66 

— 

70 

IT14 

Col. 

71 

- 

75 

IT15 

Format  ' 

1615 : 

Col. 

1 

— 

5 

NS 

Col. 

6 

- 

10 

NBM 

Col. 

11 

- 

15 

NMT 

Col. 

16 

- 

20 

NTHC 

Col. 

21 

- 

25 

NPSID 

Col. 

26 

.  - 

30 

NCASE 

CARD  04  Format  8F10.0: 

(1)  Pylon  length 

(2)  Speed  of  sound 

(3)  Air  density 

(4)  Pylon  conversion  angle 

(5)  Structural  damping  ratio 
of  rotor-blades 

(6)  Rotor  rpm 

(7)  Blade  fundamental  torsional 
natural  frequency 

(8)  Gravitational  acceleration 


( in) 
( ft/sec) 
(slug/ft3) 
(deg) 


(rpm) 
( rad/ sec  ) 

( in/sec^ ) 


CARD  05  Format  8F10.0: 


(1)  Blade  built-in  precone  angle  (deg) 

(2)  Blade  pitch-flap  coupling,  (£3) 

(3)  Blade  collective  pitch  at  initial  rpm  (deg) 

(4)  Blade  lateral  cyclic  pitch  (deg) 

(5)  Blade  longitudinal  cyclic  pitch  (deg) 

(6)  Aircraft  true  airspeed  (kt) 

(7)  Real  time  limit  (sec) 

(8)  Gravitational  acceleration  (in/sec2) 
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CARD  06 


CARD  07 


CARD  08 


Format  8F10.0: 

(1)  Integration  interval 

(  sec  ) 

(2)  Initial  azimuth  position  ot 
master  blade 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

Format  8F10.0: 

( deg  ) 

( 1 )  RPM  1 

( rpm  ) 

( 2  )  RPM  2 

(rpm) 

( 3  )  RPM  3 

( rpm) 

(4)  RPM  4 

(rpm) 

(5)  Longitudinal  hub  restraint  spring 
rate  for  rpm  1 

(  f t-lb/deg ) 

(6)  Longitudinal  hub  restraint  spring 
rate  for  rpm  2 

(ft- lb/deg) 

(7)  Longitudinal  hub  restraint  spring 
rate  for  rpm  3 

( ft-lb/deg ) 

(8)  Longitudinal  hub  restraint  spring 
rate  for  rpm  4 

Format  8F10.0: 

( ft-lb/deg ) 

(1)  Lateral  hub  restraint  spring 
rate  for  rpm  1 

( ft-lb/deg) 

(2)  Lateral  hub  restraint  spring 
rate  for  rpm  2 

( ft-lb/deg) 

(3)  Lateral  hub  restraint  spring 
rate  for  rpm  3 

( ft-lb/deg) 

(4)  Lateral  hub  restraint  spring 
rate  for  rpm  4 

(ft-lb/deg) 

(5)  Flapping  stop  spring  rate 

( ft-lb/deg) 

(6)  Blade  flapping  amplitude  at 
which  flapping  stop  is 
encountered 

(7) 

(8) 

(deg ) 

CARD  09  Format  8F10.0: 

(1)  Blade  collective  pitch  at  which  (deg) 

feathering  or  unfeathering  is  to  end 

(2)  Real  time  at  which  feathering  or  (sec) 

unfeathering  starts 

(3)  Real  time  at  which  feathering  or  (sec) 

unfeathering  ends 

(4) 
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(5) 

(6) 

(7) 

(8) 

CARD  10  Format  8F10.0: 


(1)  Amplitude  of  vertical  gust  (in/sec) 

(2)  Amplitude  of  head-on  gust  (in/sec) 

(3)  Frequency  of  sinusoidal  gust  (cps) 

(4)  Real  time  at  which  sharp-edged  (sec) 

gust  starts 

(5)  Real  time  at  which  sharp-edged  (sec) 

gust  ends 

(6)  Real  time  at  which  1-cosine  shaped  (sec) 

gust  starts 

(7)  Real  time  at  which  1-cosine  shaped  (sec) 

gust  ends 

(8) 


CARD  11  Format  8F10.0: 


(1)  Wing  angle  of  attack  (deg) 

(2)  Mean  aerodynamic  chord  of  wing  (in) 

(3)  Two-dimensional  lift  curve  slope  of  wing  (/rad) 

(4)  Wing  thickness  ratio 

(5)  Wing  sweep  angle  (deg) 

(6)  Distance  between  pylon  conversion  (in) 

hinge  and  leading  edge  of  wing  at 

wing  tip 

(7)  Distance  between  shaft  hinge  and  (in) 

wing  tip  section 

(8) 


CARD  12  Format  8F10.0: 


(1)  Initial  frequency  of  cosine  function  (cps) 

for  swashplate  lateral  excitation 

(2)  Initial  frequency  of  sine  function  for  (cps) 
swashplate  lateral  excitation 

(3)  Initial  frequency  of  cosine  function  (cps) 

for  swashplate  longitudinal  excitation 

(4)  Initial  frequency  of  sine  function  for  (cps) 
swashplate  longitudinal  excitation 

(5)  Amplitude  of  cosine  function  for  (deg) 

function  for  swashplate  lateral 

excitation 

(6)  Amplitude  of  sine  function  for  swash-  (deg) 

plate  lateral  excitation 

(7)  Amplitude  of  cosine  function  for  (deg) 

swashplate  longitudinal  excitation 

(8)  Amplitude  of  sine  function  for  swash-  (deg) 

plate  longitudinal  excitation 
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CARD  13 


Format  8F10.0: 


(1) 

Rate-of-change  of  frequency 
of  cosine  function  for  swash- 
plate  lateral  excitation 

(cyc/sec^ ) 

2 

(2) 

Rate-of-change  of  frequency 
of  sine  function  for  swash- 
plate  lateral  excitation 

(cyc/sec  ) 

2 

(3) 

Rate-of-change  of  frequency  of 
cosine  function  for  swashplate 
longitudinal  excitation 

( cyc/sec  ) 

2 

(4) 

Rate-of-change  of  frequency  of 
sine  function  for  swashplate 
longitudinal  excitation 

( eye/ sec  ) 

(5) 

Real  time  at  which  both  excita¬ 
tions  of  swashplate  start 

( sec  ) 

(6) 

(7) 

(8) 


CARD  14  Format  8F10.0: 

(1)  Component  in  a  cartesian 
coordinate  to  define  the  location 
of  an  aerodynamic  vane 

(2)  Component  in  a  cartesian 
coordinate  to  define  the  location 
of  an  aerodynamic  vane 

(3)  Component  in  a  cartesian 
coordinate  to  define  the  location 
of  an  aerodynamic  vane 

(4)  Span  of  aerodynamic  vane 

(5)  Chord  of  aerodynamic  vane 

(6)  Initial  frequency  for  aerodynamic 
vane  excitation 

(7)  Steady  angle-of-attack  of  aero¬ 
dynamic  vane,  relative  to  pylon 

(8)  Amplitude  of  angle-of-attack  oscilla¬ 
tion  of  aerodynamic  vane,  relative 

to  ,  /Ion 


( in) 

(in) 

(in) 

(in) 

(in) 

(cps) 

(deg) 

(deg) 


CARD  15  Format  8F10.0: 

(1)  Real  time  at  which  excitation  of 
aerodynamic  vane  starts 

(2)  TVo-dimensional  lift  curve  slope 
for  aerodynamic  vane 

(3)  Rate-of-change  of  the  excitation 
frequency  for  aerodynamic  vane 

(4) 

(5) 

(6) 

(7) 

(8) 


( sec  ) 
(/deg) 
( cyc/sec^  ) 
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CARD  Lb* 

CARD  17* 

CARD  18* 

CARD  19* 

CARD  20* 

CARD  21* 

CARD  22* 

CARD  23* 

CARD  24* 

CARD  25* 

CARD  26 

CARD  27 


Format  8F10.0: 

Blade  build-in  twists  (deg) 

Format  8F10.0: 

Radial  distances  to  blade  discrete  masses  (in) 

Format  8F10.0: 

Discrete  blade  masses  (chug)** 

Format  8F10.0: 

Blade  mass  unbalances  (in) 

Format  8F10.0: 

Blade  torsional  mass  moments-of-inertia  (chug-in^) 


Format  8F10.0: 

Distances  from  blade  feathering  axis  (in) 

to  blade  aerodynamic  center 

Format  8F10.0: 

Distances  from  blade  aerodynamic  center  (in) 

to  blade  center  of  gravity 

Format  8F10.0: 

Distances  from  blade  aerodynamic  center  (in) 

to  blade  elastic  axis 

Format  8F10.0: 

Blade  chord  lengths  (in) 

Format  8F10.0: 

Blade  aerodynamic  segment  lengths  (in) 

Format  13F6.0,  2X: 


Blade  root  collective  pitches  (deg  &  rad/sec) 

followed  by  rotor  angular 

velocities 

Col.  1  -  80,  alphameric. 

Title  card  to  identify  blade 
first  normal  mode 


CARD  28*  Format  13F6.0,  2X: 

Blade  natural  frequencies  for  (rad/sec) 

first  normal  mode 

*  Indicates  more  than  one  card  may  be  required  for  these  data. 

**A  unit  of  mass.  The  corresponding  unit  of  acceleration  is 
in  inches  per  second2. 
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CARD  29 


CARD  30* 

CARD  31* 

CARD  32 
CARD  36 
CARD  37 
CARD  41 
CARD  42 
CARD  46 
CARD  47 

CARD  48 

CARD  49 


*Indicates 


Format  1615: 


Col. 

1  - 

5 

LC  ( 1 ) 

Col. 

6  - 

10 

LC  ( 2 ) 

Col. 

11  - 

15 

LC  ( 3  ) 

Col. 

16  - 

20 

KDR 

Format  13F6.0,  2X: 

Normalized  blade  flexible 
out-of-plane,  in-plane,  and 
torsional  modes 

Format  9F8.0,  8X: 

Blade  normalized  bcamwise  bending 
moments  and  chordwisc  bending  moments 

Same  as  cards  27-31,  except  now  for 
the  second  normal  mode 


Same  as  cards  27-31,  except  now  for 
the  third  normal  mode 


Same  as  cards  27-31,  except  now  for 
the  fourth  normal  mode 


( in/in , 
in/in , 
rad/i n ) 


( in-lb/in) 


Format  13F6.0,  2X: 

Initial  conditions  of  normal 
coordinates  for  flexible  blades 

Format  1615: 

Col.  1-5  NSW 

Col.  6-10  NSHP 

Col.  11  -  15  NPDD 

Format  8F10.0: 


( chug  ) 
(in) 


(in) 


( in) 


more  than  one  card  may  be  required  for  these  data 


I 

-t~i  ■  rmunna  ~T>  iWiiTIMliMWI 


(1)  Wing  structural  damping  ratio 

(2)  Pylon  mass 

(3)  Component  of  distance  from  pylon 
conversion  axis  to  pylon  center  of 
gravity  in  cartesian  coordinates 

(4)  Component  of  distance  from  pylon 
conversion  axis  to  pylon  center  of 
gravity  in  cartesian  coordinates 

(5)  Component  of  distance  from  pylon 
conversion  axis  to  pylon  center  of 
gravity  in  cartesian  coordinates 


CARD  50 

CARD  51* 

CARD  52* 

CARD  53* 

CARD  54* 

CARD  55* 

CARD  56* 

CARD  57* 

*Indicates 


(6) 

Hub  mass 

(chug ) 

(7) 

Mean  distance  from  wing  aero¬ 
dynamic  center  to  its  elastic 

(in) 

axis 

r\ 

(8) 

Pylon  pitching  mass  moment-of- 
inertia  about  its  center  of 

( chug-i rr) 

gravity 

Format  8FI0.0: 

Initial  conditions  for  wing  normal 
coordinates 

Format  8F10.0: 

Wing  chordwise  normal  modes 

Format  8F10.0: 

Wing  beamwise  normal  modes 

Format  8F10.0: 

Wing  torsional  normal  modes 

Format  8F10.0: 

Discrete  wing  masses, 

wing  torsional  mass  moments-of-inertia , 

wing  mass  unbalances,  and 

wing  aerodynamic  segment  lengths 

Format  8F10.0: 

Wing/pylon  natural 
frequencies ; 

normalized  F/A  deflections, 
lateral  deflections,  and 
vertical  deflections  at  CA 


( in/ in) 


( in/in) 


( rad/in ) 


(chug, 
chug-in2 , 
in, 
in) 


(rad/sec , 
in/in , 
in/ in , 
in/ in , 


Normalized  pylon  yawings, 
pitchings,  and 
effective  masses 

Format  5D15.0: 

Elements  of  stiffness  matrix  at  a 
specified  wing  buttline 

Format  8F10.0: 

Locations  of  pylon  accelerometers, 
relative  to  pylon  conversion  axis 


rad/ in , 
rad/ in , 
chug  ) 


(  in) 


more  than  one  card  may  be  required  for  these  data 
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2 .  Guide  to  Input  Format 

CARD  OL  This  card  is  used  to  identify  the  output. 

CARD  02  This  card  is  used  to  specify  the  options  desired. 

The.  options  selected  will  depend  on  the  configura¬ 
tion  simulated  and  on  the  users  objectives.  It  is 
pointed  out  tliat  much  of  the  other  data  which  must 
be.  supplied  is  dependent  to  some  extent  on  the 
options  selected.  Any  discrepancy  between  the  data 
requirements  of  the  options  selected  and  the  data 
supplied  will  cause  the  program  to  terminate. 

IT1  allows  the  hub  restraint  springs  used  with 
a  gimbaled  rotor  to  be  simulated. 

If  IT1  =  0,  hub  restraint  springs  are 
not  simulated. 

If  IT1  =  1,  hub  restraint  springs  are 
simulated. 

IT2  is  used  to  feather  or  unfeather  the  blades 
to  start  or  stop  the  proprotor. 

If  IT2  =  0,  the  blades  will  be  feathered 
or  unfeathered. 

If  IT2  =  1,  the  blade  collective  pitch 
will  be  constant. 

IT3  is  used  to  simulate  a  discrete  gust 

encounter  or  an  oscillation  of  a  wind 
tunnel  airstream  such  as  in  the  NASA- 
Langley  16-foot  TDT. 

If  IT3  =  0,  no  gust  will  be  encountered. 

If  IT3  =  1,  a  sharp  edged  gust  will  be 
encountered . 

If  IT3  =  2,  a  1-cosine  gust  of  specified 
length  will  be  encountered. 

If  IT3  =  3,  a  continuous,  sinusoidal  gust 
is  applied. 

IT4  allows  the  program  to  be  run  with  or 
without  the  blade  elastic  degrees  of 
freedom . 
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It  IT4  =  0,  the  blade-  clastic  degrees  of 
freedom  will  be  simulated  (cards  26-47  must 
be  supplied). 

It  IT4  =  1,  the  blades  are  rigid  (cards 
26-47  must  be  deleted). 

IT5  allows  the  program  to  be  run  with  or 

without  the  wing/pylon  elastic  degrees 
of  freedom. 

If  IT5  =  0,  the  wing  degrees  of  freedom 
are  deleted  (cards  48  through  end  of  deck 
must  be  deleted). 

If  IT5  =  1,  the  wing/pylon  degrees  of 
freedom  are  simulated  (cards  48  through 
end  of  deck  must  be  supplied). 

IT6  gives  the  user  a  choice  as  to  machine  plotted 
data . 

If  IT6  =  0,  plots  will  not  be  generated. 

If  IT6  =  1,  time,  histories  will  be 
plotted  for  rotor  rpm,  longitudinal  and 
lateral  flapping,  rotor  H-force,  rotor  Y- 
force,  rotor  thrust,  and  torque. 

If  IT6  =  2,  time  histories  will  be  plotted 
for  the  inplane  and  out-of-plane  deflections 
of  the  master  blade,  beamwise  and  chordwise 
bending  moments  at  the  first  specified  blade 
station,  beamwise  and  chordwise  moments  at 
the  second  specified  blade  station  and 
flapping  in  the  rotating  system. 

i 

If  IT6  =  3,  time  histories  will  be  plotted 
for  the  wing  tip  beamwise,  chordwise  and, 
torsional  deflections,  vertical  accelerations 
of  the  first  two  specified  accelerometers, 
and  the  axial  and  lateral  accelerations 
at  the  first  accelerometer. 

IT7  allows  the  aerodynamic  interference  between 
the  wing  and  rotor  to  be  simulated. 

If  IT7  =  0,  the  aerodynamic  interference 
will  be  zero. 

If  IT7  =  1,  the  aerodynamic  interference 
will  be  simulated. 
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ITS 


allows  the  user  to  excite  the  system  using 
either  the  swashplate  or  an  aerodynamic  vane 
Located  on  the.  pylon. 

It  ITS  =  0,  there  is  no  excitation. 

It  ITS  =  I,  the  swashplate  will  be  excited 
cyclically. 

If  IT8  =  2,  the  aerodynamic  vane  will  be 
osciLlated  in  pitch. 

IT9  allows  the  user  to  collect  peak  values  for 
several  variables  throughout  the  time 
history.  This  option  prints  out  at  the 
end  of  the  output  data  the  maximum  and 
minimum  values  of  master  blade  flapping 
(rotating  system),  beamwise  and  chordwise 
deflections  of  the  wing  tip,  and  the  first 
normal  coordinate  of  the  wing/pylon  system. 

If  IT9  =  0,  peaks  will  not  be  printed  out. 

If  IT9  =  1,  peaks  are  printed  out. 

IT10  allows  the  accelerations  at  specified 

pylon  stations  to  be  calculated  and  printed 
out . 

If  IT10  =  0,  the  accelerations  will  not  be 
printed. 

If  IT10  =  1,  the  accelerations  at  up  to 
nine  specified  locations  will  be  printed. 

If  IT5  =  0,  this  option  is  ignored  by  the 
program. 

IT11  allows  wing  bending  moments  to  be  calculated. 

If  IT11  =  0,  wing  bending  moments  will  not 
be  calculated. 

If  IT11  =  1,  wing  bending  moments  at  a 
specified  station  will  be  calculated  and 
printed  out. 

IT12  allows  the  user  a  choice  of  rigid-body 
freedoms  of  the  proprotor. 

If  IT12  =  0,  a  gimbaled  semirigid  rotor, 
with  rotational  freedom  (variable  rpm) 
is  simulated. 
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It  ITL2  =  I,  the  blade  rigid-body  flapping 
is  locked  out  to  allow  simulating  ot  hinge¬ 
less  or  articulated  rotors.  Rotor  rprri  is 
a  variable. 


If  IT12  =  2,  the  rotor  rpm  is  constant 
(flapping  is  free). 

If  IT12  =  3,  the  rigid-body  flapping  and 
rotational  freedoms  are  locked  out. 

IT13  allows  the  effective  mass  of  the  wing/pylon 
modes  to  be  calculated  by  the  program. 

If  IT13  =  0,  the  wing/pylon  normal  modes 
effective  masses  are  calculated  internally. 

If  IT13  =  1,  input  values  of  the  wing/pylon 
effective  masses  are  used  (card  55). 

IT14  is  used  to  overcome  a  numerical  problem  often 
encountered  with  time  history  analyses.  In 
some  cases  pylon  accelerations  with  fre¬ 
quencies  the  same  as  one-half  the  reciprocal 
of  the  integration  interval  have  been 
encountered.  By  averaging  the  accelerations 
to  be  calculated  these  high  frequency 
accelerations  can  be  eliminated. 

If  IT14  =  0f  no  averaging  takes  place. 

If  IT14  =  1,  calculated  accelerations  are 
averaged. 

IT15  is  not  used  at  present. 

CARD  03  is  used  to  read  in  parameters  which  specify  the 
amount  of  the  input  data  to  follow. 

NS  is  the  number  of  blade  lumped  mass  stations. 

0  <  NS  <  10 

NBM  is  the  number  of  normal  blade  modes  used. 

0  <  NBM  <  4.  NBM  is  used  only  if  IT4  =  0. 

NMT  is  the  number  of  blade  stations  at  which 

blade  bending  moments  are  to  be  calculated. 

0  <  NMT  <  4.  NMT  is  used  only  if  IT4  =  0. 

NTHC  is  the  number  of  blade  collective  pitch 
angles  used  in  the  blade  normal  mode 
tables.  0  <  NTHC  <  6 
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NPSID  is  the  number  of  rotor  rpm  used  in  the 

blade  normal  mode  tables.  0  <  NPSID  ^  6 

NGASE  number  of  cases  in  a  single  computer  run. 

It  is  read  in  from  the  data  deck  of  the 
first  case  only.  If  NGASE  is  left  blank, 
the  program  will  assume  there  is  only  one 
case. 

CARD  04  contains  descriptive  data. 

ZH  is  the  pylon  length;  the  distance  from  the 

pylon  conversion  axis  to  the  rotor  hub. 

S  is  the  speed  of  sound;  must  be  greater  than 

zero. 

RHO  is  the  air  density.  If  it  is  input  as  zero, 
computation  of  rotor  aerodynamics,  wing 
aerodynamics,  pylon  aerodynamic  excitation, 
and  the  wing/rotor  aerodynamic  interference 
is  not  made.  If  user  does  not  wish  to 
include  rotor  aerodynamic  forces,  zeros 
should  be  input  on  card  25.  If  user  does 
not  wish  to  include  wing  aerodynamics, 
zeros  should  be  input  for  the  wing  aero¬ 
dynamic  segment  lengths  on  card  54. 

TPG  is  the  pylon  conversion  angle,  measured 
relative  to  flight  path.  It  is  zero  for 
a  proprotor  in  propeller  mode,  90  degrees 
for  helicopter  mode. 

BZETA  is  the  blade  structural  damping  ratio, 
ignored  by  the  program  if  IT4  =  1. 

RPM  is  initial  rotor  rpm.  rpm  £  0 

OMT  is  blade  fundamental  torsional  natural 

frequency.  OMT  >  0.  OMT  is  used  in  the 
calculation  of  blade  pitch  damping 
coefficients . 

G  is  the  gravitational  acceleration  normal 

to  the  wing.  G  >  0 

CARD  05  ag  is  the  rotor  precone  angle. 

6_  is  the  blade  pitch-flap  coupling.  Positive 

if  blade  pitch  decreases  when  the  blade 
flaps  upward. 
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is  the  blade  initial  collective  pitch. 
Corresponds  to  initial  rotor  rpm. 

is  the.  blade  lateral  cyclic  pitch  (sign 
convention  is  that  given  in  Reference  1). 

is  the  blade  longitudinal  cyclic  pitch. 

U  is  the  velocity  along  flight  path.  Must 

be  greater  than  or  equal  to  zero. 

REALTM  the  real  time  required  for  the  case. 

Program  execution  is  contained  in  the  time 
history  loop  until  REALTM  is  reached.  To 
avoid  waste  of  machine  time  due  to  input 
errors  the  program  is  terminated  if  the 
rotor  rotational  speed  is  below  -2  rad/sec. 

GY  Gravitational  acceleration  parallel  to  the 
wing  span.  For  example,  GY  is  used  for 
a  wind  tunnel  floor-mounted  semispan 
model  wing. 

CARD  06  DTMIN  is  the  time  increment  used  in  the  time 

history  calculation.  The  user  should  not 
attempt  to  use  a  DTMIN  which  is  larger 
than  the  equivalent  of  /■.  =  15  degrees. 

A  value  of  DTMIN  =  0.005  second  is  picked 
by  the  program  if  it  is  not  specified. 
Experience  suggests  the  use  DTMIN  =  0.005 
second  for  runs  where  the  blade  normal 
modes  are  included.  Otherwise  the  following 
DTMINs  may  be  used: 

For  rotor  rpm  £  300  DTMIN  =  0.006  second 

300  rotor  rpm  £  200  DTMIN  =  0.007  second 

200  rotor  rpm  £  100  DTMIN  =  0.01  second 

rotor  rpm  <  100  DTMIN  =  0.012  second 

Since  averaging  process  is  taken  if  IT14  =  1, 
integration  interval  will  then  be  twice  of 
DTMIN.  This  doubled  value  should  not  be 
larger  than  the  equivalent  of  A'.  =  15  degrees. 

PSIE  is  the  initial  azimuth  position  of  the 
master  blade. 
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CARD  07 


The  data  on  this  card  can  be  used  to  simulate  a 
nonisotropic,  variable  rate,  hub  restraint  spring. 
The  rpm  range  is  divided  into  three  parts,  and  the 
variation  of  spring  rate  is  assumed  to  be  linear 
over  these  three  regions,  see  Figure.  2. 


RPM 


Figure  2.  Hub  Restraint  Spring  Rate-rpm  Relationship. 


The  first  four  floating 
RPM3,  and  RPM4  give  the 
rate  is  specified.  The 
AKi3,  and  AK14 ,  specify 
rate  at  these  rpm. 


point  numbers ,  RPM1 ,  RPM2 , 
rpm  at  which  the  spring 
last  four,  AKu  ,  AK12, 
the  longitudinal  hub  spring 


CARD  08  This  card  is  used  in  conjunction  with  Card  07.  The 
first  four  floating  point  numbers  BK11,  BK12 ,  BK13, 
and  BK14  give  the  lateral  hub  spring  rates  at  the 
four  rpm  specified  on  Card  07.  The  fifth,  STOPK, 
is  a  hub  spring  rate  simulating  flapping  stops.  The 
sixth,  STOPTH,  is  the  blade  flapping  angle  at  which 
the  blade  (or  hub)  contacts  the  flapping  stop.  If 
the  longitudinal  hub  spring  rate  on  Card  07  is 
constant,  and  the  lateral  hub  spring  rate  on  Card 
08  is  also  constant;  it  is  then  suggested  that 
zeros  be  input  for  each  rpm  to  avoid  an  inter¬ 
polation  requirement. 
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CARD  09  This  card  controls  blade  feathering.  It  is  ignored 
if  IT2  =  1. 

THCMAX  Tti is  is  the  blade  collective  pitch  at 

which  feathering  or  unfeathering  stops. 

TMST  Real  time  at  which  feathering  or 

unfeathering  starts . 

TMED  Real  time  at  which  feathering  or 

unfeathering  ends. 

CARD  10  This  card  specifies  the  magnitude  and  time  a  gust 
is  encountered.  The  magnitude  and  time  relation¬ 
ships  arc  illustrated  in  Figure  3. 

(a )  Step  Gust 


Figure  3.  Gust  Magnitude  and  Time  Relationships. 
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Oas/Ni)  (oas/Ni)  isn0 


(b) 


L  -  cosine  gust 


V 

gus 


t 


=  a 


(c)  Sinusoidal  gust 


Figure  3.  Concluded. 
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CARD  LL 


This  card  is  used  to  specify  wing  aerodynamic 

and  wing/rotor  aerodynamic  interference,  parameters. 

If  IT5  =  0  and  IT 7  =  0,  the  data  on  this  card 

will  not  be  used. 

ALW  is  the  wing  angle  of  attack.  Positive 

wing  Leading-edge  up. 

WC  is  the  wing  mean  aerodynamic  chord,  must 

be  greater  than  zero. 

SLC  is  the  two-dimensional  lift  curve  slope 

for  wing.  Stall  is  not  considered. 

DEL  is  the  ratio  of  wing  thickness  to  wing 

mean  chord. 

ALAM  is  the  wing  sweep  angle,  positive  forward. 

AL  Illustrated  in  Figure  3. 

YCA  Illustrated  in  Figure  3. 


PROPROTOR 


Figure  4.  Rotor/Wing  Geometrical  Relations. 
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CARD  12 


This  card  and  Card  13  arc  used  lor  swashplate 
cyclic  excitation.  For  ITS  =  1,  the  position  oi 
the  swashplate  is  given  by: 

tor  t  -  t 


AL(t)  =  Al 
B  L ( t )  = 

for  tit 


AL(t)  =  Al 

+  A. 
lc 

cos 

['acf'1  1 

+  Als 

sin 

as 

|t-t  )' 

BL(t)  =  Bl 

+  B. 

l.c 

cos 

.-bcf-1  ). 

+  Als 

sin 

bs 

l-1  )] 

where  t  is  the  time  excitation  initiated. 
And , 


jj 


ac 


u) 

as 


be 


CO 


bs 


J  acO 

asO 

"bcO 

^bsO 


ac 


at 


+ 


dc 

as 

at 


M:) 


The  eight  floating 
input  in  the  order 

Alc’  Als’  Blc ’  and 


point  numbers  on  this  card  are 
that:  o)ac0,  (oag0,  o)bc0,  ^bsQ, 

Bls* 


CARD  13  Inputs  are  in  the  following  order:  /at, 

a^as/at,  a^bc/at,  aobg/at,  and  t„.  a 

CARD  14  The  data  on  this  card  and  Card  15  are  used  for  an 
aerodynamic  vane  mounted  on  the  pylon  to  excite  the 
wing/pylon  normal  modes.  The  location  and  dimen¬ 
sions  of  the  exciter  can  be  specified.  The  angle 
of  attack  of  the  vane  is  calculated  by 
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e  ( 1  "l  e ) 

in  which  (  pc  is  the  pylon  conversion  any  1  e ,  Qe  Jnd 
oqs  are  the  steady  and  oscillatory  angle  ol  attack 
ot  the  vane  relative  to  the  pylon  respectively,  tc 
is  the  time  when  excitat  ion  ot  the  vane  begins, 
and  t.  is  the  exciting  frequency  given  by: 


a  = 
e 


pc 


+  t. 


'Qe 


°0s  sln 


o 

e  eO  i)  t  ( *  *  e  ) 


UE  \ 

a  re 

the 

Car 

tesian 

coordinates  used  to 

loc 

ate 

the 

aerodyn 

lamic  vane  in  the  pylon 

VE  ) 

coo 

rdin 

ate 

system. 

The  intersection  of 

the 

qua 

rter 

chord 

and  midspan  ol  the  vane 

WE  f 

is 

used 

as 

the  reference  point. 

SPNE 

is 

the 

span 

of  the 

vane . 

CDE 

is 

the 

aerodynamic 

chord  of  the  vane. 

OME 

is 

the 

eO 

in  the 

expression  for 

e 

ALX 

is 

the 

Oe 

in  the 

equation  for 

Pos 

itiv 

^  if 

the  leading  edge  of 

the 

vane  is 

up . 

ALS 

is 

the 

"Os 

in  the 

equation  for  cc. 

It  is  positive  if  leading  edge 
of  the  vane  is  up. 


CARD  15  This  card  is  a  continuation  of  Card  14. 


TEX 

is 

the 

time 

excitation  begins 

( te ) . 

CLXO 

is 

t  he 

lift 

curve  slope  of  the 

vane . 

OMES 

is 

du 

e 

the 

/at. 

rate 

of  change  of  vane 

frequency , 

CARD  16*  Card  16  through  25  contain  blade  geometry  and 
inertial  data. 

THTtI)  is  the  blade  built-in  twist  for  a  given 

blade  segment.  THT  is  read  in  from  blade 
root  to  blade  tip.  Washout  is  negative 
thrust  (I  =  1  to  NS). 

CARD  17*  R(I)  is  the  radial  distance  for  the  blade- 

lumped  masses.  R(I)  is  read  in  from  root 
to  tip  (I  =  1  to  NS ) . 

-'Indicates  more  than  one  card  may  Jdc.  required  for  these  data. 
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CARD  L8*  AM( I )  is  the  blade  Lumped  masses.  AM( I  )  is  read 

in  from  root  to  tip  (I  =  1  to  NS). 

CARD  19*  YB(I)  is  the  blade  mass  unbalance.  YB(I)  is 

read  in  from  root  to  tip  (I  =  1  to  NS). 

CARD  20*  AIC(I)  is  the  blade  torsional  mass  moment  of 

inertia.  AIC(I)  is  read  in  from  root 
to  tip  (I  =  1  to  NS). 

CARD  21  EPF(I)  is  the  distance  from  the  blade  feathering 

axis  to  the  blade  aerodynamic  center. 
Positive  if  AC  is  ahead  of  the  feathering 
axis.  Input  is  read  in  from  root  to  blade 
tip  (I  =  1  to  NS). 

CARD  22*  EPC(I)  is  the  distance  from  the  blade  aerodynamic 

center  to  the  blade  center  of  gravity. 
Positive  if  the  center  of  gravity  is  ahead 
of  the  aerodynamic  center  (I  =  1  to  NS). 

Read  in  from  blade  root  to  blade  tip. 

CARD  23*  EPE(I)  is  the  distance  from  the  blade  aerodynamic 

center  to  the  blade  elastic  axis.  Positive 
if  aerodynamic  center  is  ahead  of  the 
elastic  axis  ( I  =  1  to  NS).  Read  in  from 
blade  root  to  blade  tip. 

CARD  24*  C(I)  is  the  blade  segment  chord  length.  Input 

from  root  to  tip  (I  =  1  to  NS). 

CARD  25*  DR(I)  is  the  width  of  the  blade  segments.  A 

zero  value  of  DR(I)  can  be  used  to  eliminate 
the  blade  aerodynamics.  By  varying  the 
segment  lengths  at  the  blade  root  and  tip, 
blade  root  cut-out  and  tip  loss  may  be 
simulated . 

CARDS  26-47  are  used  to  input  the  blade  normal  modes.  The 
information  on  Card  26  is  used  for  all  of  the 
modes.  Cards  27  through  31  contain  data  describing 
the  first  blade  mode.  Cards  32  through  36,  37 
through  41,  and  42  through  46  contain  data  for  the 
second,  third  and  fourth  modes  respectively.  If 
IT4  =  1,  cards  26  through  47  should  be  deleted  from 
the  data  deck.  If  IT4  =  0,  then  NBM  on  Card  03  will 
determine  how  many  data  cards  must  be  deleted  from 
the  data  deck. 

'’"Indicates  more  than  one  card  may  be  required  for  these  data. 


The  blade  normal  mode  data  are  supplied  in  a  three-dimensional 
array  format  as  Illustrated  in  Figure  5. 


Figure  5.  Three-Dimensional  Array  Format  for 
Blade  Normal  Modes  Data. 

CARD  26  contains  the  collective  pitch  and  rotor  angular 

velocity  arguments.  The  number  of  collective  pitch 
values  and  rotor  angular  velocity  values  to  be 
read  is  specified  by  NTHC  and  NPSID  on  Card  03, 
respectively.  Should  the  collective  pitch  or  rpm 
fall  outside  the  tabulated  domain,  extrapolation 
will  not  be  made.  Instead,  the  program  will  use 
the  closest  value  obtainable  from  the  table,  i.e. 
a  value  on  the  boundary  of  the  tabulated  domain. 

CARD  27  is  used  to  identify  the  first  input  blade  mode, 

i.e.  collective  or  cyclic,  flatwise  or  edgewise  etc. 

CARD  28*  contains  the  natural  frequencies,  in  radians  per 

second,  for  the  first  normal  blade  mode.  The  fre¬ 
quencies  are  read  in  first  as  a  function  of  the 
collective  (pitch  for  the  lowest  rpm),  then  for 
successively  higher  rpm. 

'^Indicates  more  than  one  card  may  be  required  for  these  data. 
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CARD  29 


i 


alLows  the  user  to  specify  if  the  shapes  of  the 
normal  modes  are  functions  of:  (1)  both  rotor 
rpm  and  collective,  (2)  collective  pitch  only,  and 
(3)  rpm  only. 

If  LC(L)  =  1,  means  the  mode  shapes  are  functions 
of  rotor  rpm  and  collective  pitch. 

If  LC(L)  =  2,  means  the  mode  shapes  are  functions 
of  collective  pitch  only. 

If  LC(L)  =  3,  means  the  mode  shapes  are  functions 
of  rpm  only. 

If  LC(L)  =  4,  means  the  mode  shapes  do  not  vary 
with  rpm  or  collective  pitch. 

There  are  two  types  of  mode  shapes  associated  with 
gimballed  rotors,  collective  modes,  and  cyclic 
modes.  The  former  are  polar  symmetric  modes;  the 
latter  polar  antisymmetric  modes.  KDR  on  Card  29 
specifies  whether  the  input  mode  is  a  collective 
or  cyclic  mode.  KDR  is  used  only  if  IT12  equals  0 
or  2 . 

If  KDR  =  0,  a  collective  mode  is  input. 

If  KDR  =  1,  a  cyclic  mode  is  input. 


CARD  30*  contains  the  matrix  of  out-of-plane,  in-plane,  and 
torsional  components  for  the  blade  normal  mode. 

The  out-of-plane  component  is  read  in  first,  followed 
by  the  inplane  and  torsional  components.  Each 
component  is  read  in  from  root  to  tip,  first  versus 
collective  pitch,  then  versus  rpm. 


CARD  31*  contains  the  matrix  of  beamwise  and  chordwise 
normalized  bending  moment  for  the  first  normal 
mode.  The  beamwise  moments  are  read  in  first, 
followed  by  the  'hr  iwise  moments.  Each  component 
is  read  in  fi^st  m  a  s  collective  pitch,  then 
versus  rpm.  ism’  0  d  03  specifies  the  number 
of  blade  stations  w.  .  a  moment  data  will  be  read 
in. 

CARDS  32*  to  36*,  37*  to  41*,  and  42*  to  46*  are  arranged 
in  the  same  order  as  Cards  27  through  31.  The 
data  on  cards  32  to  36  are  for  the  second  blade 
mode  and  should  be  deleted  if  NBM  on  Card  03  <  1. 

-'Indicates  more  than  one  card  may  be  required  for  these  data. 
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The  data  on  cards  37  to  41  are  for  the  third  blade 
mode  and  should  be  deleted  if  NBM  on  Card  03  ^  2. 

The  data  on  cards  42  to  46  are  for  the  fourth  blade 
mode  and  are  deleted  if  NBM  on  Card  03  £  3. 

CARD  47  contains  the  initial  conditions  for  the  normal 

blade  modes.  The  initial  conditions  for  the  first 
mode  of  each  blade  are  read  in  first,  followed  by 
those  for  the  second  blade  mode,  etc. 

CARD  48  is  used  to  specify  parameters  which  indicate  the 
amount  of  the  input  data  to  follow. 

NSW  is  the  number  of  wing  lumped  mass  stations. 

0  <  NSW  <  10.  NSW  is  read  in  only  if 
IT5  =  1. 

NSHP  is  the  number  of  normal  modes  used  for  the 
wing/pylon  system.  0  <  NSHP  t  4.  NSHP  is 
read  in  only  if  IT5  =  1. 

NPDD  is  the  number  of  pylon  accelerometers. 

0  <  NPDD  <  9.  NPDD  is  used  only  if 
IT5  =  1  and  IT10  =  1. 

CARD  49  contains  data  for  the  wing/pylon  system. 

WZETA  is  the  wing  structural  damping  ratio.  A 
zero  value  is  acceptable. 

PM  is  the  lumped  pylon  mass  (complete  nacelle 

mass  less  the  rotor  group)  PM  >  0. 

2P  is  the  location  (Z  direction)  of  the  pylon 

center  of  gravity  along  the  shaft  center- 
line  from  the  conversion  axis  (see  Figure  6). 

YP  is  the  lateral  location  (Y  direction)  of 

the  pylon  center  of  gravity  with  respect  to 
the  shaft  axis  (see  Figure  6). 

XP  is  the  vertical  location  (X  direction)  of 

the  pylon  center  of  gravity  with  respect  to 
the  shaft  axis  (see  Figure  6). 

HM  Concentrated  hub  mass  (rotor  group  mass  less 

the  blade's  mass)  HM  >  0  (see  Figure  6). 

WEPE  is  the  mean  distance  from  the  wing  aero¬ 
dynamic  center  to  the  wing  elastic  axis. 
Positive  if  the  aerodynamic  center  is  ahead 
of  the  elastic  axis. 
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AND  INERTIA 


Figure  6.  Pylon  Geometries. 


WIPP  is  the  pylon  pitching  mass  moment-of- 

inertia  about  the  pylon  center  of  gravity. 

CARD  50  contains  SHPE(J) ,  the  initial  conditions  for  the 
wing  normal  modes  (J  =  1  to  NSHP) . 

CARD  51*  contains  the  chordwise  components  of  the  wing  normal 
modes.  They  are  read  in  from  wing  root  to  wing  tip 
for  first  mode,  then  for  the  second  mode,  etc.  The 
number  of  the  discrete  wing  stations  is  specified 
by  NSW  on  Card  48. 

CARD  52*  contains  the  beamwise  components  of  the  wing 

normal  modes.  They  are  read  in  the  same  order  as 
the  chordwise  components. 


♦Indicates  more  than  one  card  may  be  required  for  these 
data . 
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CARD  53* 


contains  the  torsional  components  ot  the  wing  normal 
modes.  They  are  read  in  the  same  order  as  the 
chordwise  components. 

CARD  54  contains  the  wing  lumped  inertial  properties.  The 
masses  from  wing  root  to  tip  are  read  in  first 
followed  by  the  torsional  mass  moment s-of- inert i a 
about  the  wing  elastic  axis,  and  the  wing  mass 
unbalances.  The  wing  segment  widths  are  read  in 
last.  The  number  of  discrete  wing  stations  for 
which  data  must  be  supplied  is  specified  by  NSW 
on  Card  48.  Wing  mass  unbalances  are  considered 
positive  if  wing  center  of  gravities  are  in  front 
of  the  wing  elastic  axis.  Wing  aerodynamic  segment 
lengths  can  be  set  to  zero  in  order  to  zero  out  wing 
aerodynamic  forces. 

CARD  55*  contains  the  natural  frequencies  and  pylon  mode 
shapes  for  the  wing  normal  modes.  The  number  of 
modes  is  specified  by  NSHP  on  Card  48.  The 
natural  frequencies  are  read  in  first,  in  succession 
for  each  mode  followed  by  the  normalized  pylon  fore 
and  aft  deflections,  the  normalized  pylon  lateral 
deflections,  the  normalized  vertical  deflections, 
the  normalized  pylon  yaw  angles,  the  normalized 
pylon  pitch  angles,  and  the  wing/pylon  effective 
masses.  The  wing/pylon  effective  masses  will  be 
calculated  by  the  program  if  IT13  =0.  In  this 
case,  the  input  values  for  the  wing/pylon  effective 
mass  are  not  used. 

CARD  56*  contains  a  modified  wing  stiffness  matrix  by  rows. 

The  wing  moments  are  calculated  using  the  relation¬ 
ship  between  force  and  normal  coordinates  given  by: 

{F}  =  [k]  {x} 

in  which  {F }  is  the  force  matrix,  {x}  the  normal 
coordinates  matrix,  and  [k]  the  modified  wing 
stiffness  matrix  obtained  by  premultiplying  the 
mode  shape  matrix  by  a  stiffness  matrix.  The  order 
of  [k]  is  determined  by  the  number  of  wing  normal 
modes,  i.e.  the  order  of  [k]  =  10  x  NSHP  on  Card  48. 


'"Indicates  more  than  one  card  may  be  required  for  these  data. 
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CARD  57*  is  used  to  specify  the  location  of  up  to  nine  (9) 

accelerometers  located  on  the  pylon.  The  locations 
are  referred  to  the  pylon  coordinate  system.  The 
X,  Y ,  Z  components  are  read  in  first,  followed 
sequentially  by  those  for  the  second  accelerometer 
etc.  The  number  of  accelerometers  is  specified  by 
NPDD  on  Card  48. 

3 .  Output  Format  and  Guide 

The  output  of  ARAP06  consists  of  two  parts:  (1)  the  input 
data  and  (2)  a  time  history  of  the  system  displacements, 
accelerations,  and  loads.  Machine  (CALCOMP)  plots  of  selected 
variables  may  be  obtained  if  desired. 

The  print-out  of  the  input  data  is  in  order  of  the  flight 
condition  being  simulated,  options,  proprotor  data,  and 
wing/pylon  data.  All  of  the  data  used  in  the  case  is  printed 
out . 

The  printed  time  history  output  may  be  conveniently  divided 
into  seven  groups.  The  first  group,  consisting  of  the  first 
two  lines,  is  always  printed.  The  printing  of  the  other 
groups  depends  upon  the  options  selected. 

a.  Proprotor  Rigid-Body  Motion  and  Forces  and  Moments 
(Rows  1  and  2J 

TIME  is  the  real  time  (sec) 

PSID  is  the  rotor  instantaneous  rotational 
speed  about  the  rotor  shaft  axis  (rpm) 

PSI  is  the  azimuth  position  of  the  master  blade, 

measured  from  X-axis  of  the  rotor-fixed 
system  (deg) 

BET-LG  is  the  longitudinal  flapping  in  rotor-fixed 
system  positive  about  negative  Y-axis  (deg) 

BET-LT  is  the  lateral  flapping  in  rotor-fixed 

system,  positive  about  positive  X-axis  (deg) 

BET(l)  is  the  total  flapping  of  the  master  blade, 
positive  if  flapping  up  (deg) 

FFX  is  the  rotor  aerodynamic  H-force  referred 

to  rotor-fixed  system.  Positive  along 
positive  X-axis  (lb) 

Indicates  more  than  one  card  may  be  required  for  these  data. 


29 


FFY 

FFZ 

WBX 

WBY 

WBZ 

FI 

FJ 

FK 

MI 

MJ 

MK 

TFFX 

TFFX 


is  the  rotor  aerodynamic  side  force 
referred  to  rotor-fixed  system.  Positive 
along  positive  Y-axis  (lb) 

is  the  rotor  aerodynamic  drag  force, 
referred  to  rotor-fixed  system.  Positive 
along  positive  Z-axis  (lb) 

is  the  X-component  of  the  gravitational 
force  on  rotor-blades,  referred  to  rotor- 
fixed  system.  Positive  along  positive 
X-axis  (lb) 


is  the  Y-component  of  the  gravitational 
force  on  rotor-blades,  referred  to  rotor- 
fixed  system.  Positive  along  positive 
Y-axis  (lb) 

is  the  Z-component  of  the  gravitational 
force  on  rotor-blades,  referred  to  rotor- 
fixed  system.  Positive  along  positive 
Z-axis  (lb) 


is  the  component  of  rotor 
in  rotor-rotating  system, 
x-axis  (lb) 

is  the  component  of  rotor 
in  rotor-rotating  system, 
y-axis  (lb) 

is  the  component  of  rotor 
in  rotor-rotating  system, 
z-axis  (lb) 

is  the  component  of  rotor 
in  rotor-rotating  system, 
x-axis  (in- lb) 

is  the  component  of  rotor 
in  rotor-rotating  system, 
y-axis  (in-lb) 

is  the  component  of  rotor 
in  rotor-rotating  system, 
z-axis  (in-lb) 


inertial 

loads 

positive 

along 

inertial 

loads 

positive 

along 

inertial 

loads 

positive 

along 

inertial 

moments 

positive 

about 

inertial 

moments 

positive 

about 

inertial 

moments 

positive 

about 

is  the  component  of  net  rotor  force 
referred  to  rotor-fixed  system.  Positive 
along  X-axis  (lb).  And 

=  FFX  +  FI  cos -  FJ  sim  +  WBX 
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TFFY  is  the  component  of  net  rotor  force 

referred  to  rotor  fixed  system.  Positive 
along  Y-axis  (lb).  And 

TFFY  =  FFY  +  FI  sin,  +  FJ  cos i  +  WBY 

TFFZ  is  the  component  of  net  rotor  force 

referred  to  rotor  fixed  system.  Positive 
along  Z-axis  (lb).  And 

TFFZ  =  FFZ  +  FK  +  WBZ 


TFMX  is  the  component  of  net  rotor  moment 

referred  to  rotor  fixed  system.  Positive 
about  X-axis  (in-lb) 


TFMY  is  the  component  of  net  rotor  moment 

referred  to  rotor  fixed  system.  Positive 
about  Y-axis  (in- lb) 


TFMZ  is  the  component  of  net  rotor  moment 

referred  to  rotor  fixed  system.  Positive 
about  Z-axis  (in-lb) 


b.  Blade  Deflections  and  Bending  Moments 


This  group  is  calculated  and  printed  only  if  IT4  =  0.  The 
bending  moments  are  written  symbolically  as  follows: 


JBMN  where  BM  stands  for  a  beamwise  moment,  J 

stands  for  the  Jth  blade,  N  stands  for  the 
Nth  blade  station  at  which  the  moment  is 
calculated  (in-lb) 

JCMN  where  CM  stands  for  chordwise  moments, 

J  and  N  have  been  defined  above  (in-lb) 

The  blade  tip  deflections  are  with  respect  to  the  preconed 
axis,  and  are  written  symbolically  as  follows: 

ZTIPJ  where  TIP  stands  for  blade  tip,  J  stands 
'  for  the  Jth  blade,  Z  indicates  a  beamwise 
deflection.  It  is  positive  if  the  blade 
is  being  bent  upward.  (in) 

XTIPJ  where  X  indicates  a  chordwise  deflection. 
It  is  positive  in  the  lag  sence,  i.e. 
compression  in  trailing  edge  (in) 

The  normal  coordinates  for  the  blade  normal  modes  are  written 
symbolically  as  follows: 
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JQRN  where  QB  stands  tor  the  normal  coordinates, 

J  indicates  the  Jth  blade,  N  stands  tor  the 
Nth  normal  mode 

c .  Deflections  and  Rotations  of  Wing/Pylon  System 
This  group  is  executed  and  printed  only  if  IT5  =  1. 

THP  is  the  pylon  instantaneous  pitch  angle, 
positive  about  Yp-axis  in  the  inertial 
fixed  system.  THP  contains  the  pylon 
steady  conversion  angle  as  well  as  the 
pitch  due  to  wing  torsional  deflections 
(rad) 

PSIP  is  the  pylon  instantaneous  yaw  angle, 

positive  about  Zp-axis  in  the  inertial 

fixed  system  (rad) 

XHB  is  the  Xp-component  of  hub  deflection, 

referred  to  reference  inertial  coordinate 
system  (in) 

XHB  =  Xo  *  zH(COS’p  cos  p'O 

in  which  X  is  Xp-component  deflection  of 
pylon  conversion  axis,  zH  is  the  total 
pylon  length  (CA  to  hub),  p  is  the  pylon 
yawing  angle ,  and  Pp  is  the  pylon  pitching 
ang  le 

YHB  is  the  Yp-component  of  hub  deflection, 

referred  to  reference  inertial  coordinate 
system  (in) 

YHB  =  Yo  *  ZH  slm:P 

ZHB  is  the  Zp-component  of  hub  deflection, 

referred  to  reference  inertial  coordinate 
system  (in) 

Z77n  =  Z  -  z„  cos'j'  sinn 
HB  o  H  p  p 

WBM  is  the  beamwise  deflection  at  the  wing-tip, 

positive  down  (in) 

WCD  is  the  chordwise  deflection  at  the  wing-tip, 

positive  forward  (in) 
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WTN 


is  t tie  torsional  deflection  at  the  wing-tip, 
positive  leading -edge  up  (rad) 


QW1 

is  the  first  normal 

coordinate  of  the 

wing/pylon  sysl 

Lem 

gw  2 

Second,  third, 

and 

fourth  normal 

gw  3 

coordinates  of 

t  he 

wing/pylon  system, 

gw4 

respectively 

d .  Wing  Loads  and  Moments 

The  execution  and  printing  of  this  group  are  controlled  by 
option  IT11.  If  IT11  =  1,  the  following  variables  are 
printed . 

FXL  is  the  chordwise  shearing  force  at  the 

left  end  of  a  wing  segment,  positive 
forward  ( lb ) 

FZL  is  the  beamwisc  shearing  force  at  the 

left  end  of  a  wing  segment,  positive 
down  (lb) 

MXL  is  the  beamwise  bending  moment  at  the 

left  end  of  a  wing  segment,  positive  if 
tension  is  produced  on  the  wing  upper 
surface  (in-lb) 

MYL  is  the  torsional  moment  at  the  left  end 

of  a  wing  segment,  positive  if  wing 
leading  edge  is  up  (in-lb) 

MZL  is  the  chordwise  bending  moment  at  the  left 

end  of  a  wing  segment,  positive  if  tension 
is  produced  on  the  wing  surface  at  the 
leading-edge  (in-lb) 

The  remaining  variables  in  this  group  are  chordwise  and 
beamwise  shearing  forces,  beamwise,  torsional,  and  chordwise 
bending  moments  for  the  right  end  of  the  wing  segment. 


e .  Pylon  Accelerations 

This  group  is  controlled  by  option  IT10.  If  IT10  =  1,  the 
computation  of  pylon  accelerations  is  executed. 


XPDD1  is  the  X-component  acceleration  of  the 
first  accelerometer,  referred  to  pylon 
coordinate  system  (in  fraction  of  g) 
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YPDDL  is  the  Y-component  acceleration  of  the 
first  accelerometer,  referred  to  pylon 

coordinate  system  (in  fraction  of  g) 

ZPDD1  is  the  Z-component  acceleration  of  the 
first  accelerometer,  referred  to  pylon 

coordinate  system  (in  fraction  of  g) 

The  remaining  variables  in  this  group  arc  component  accelera¬ 
tions  for  the  accelerometers,  --NPDD. 

f .  Peak  Values  of  Selected  Variables 

The  printing-out  of  this  group  is  controlled  by  option  IT9. 

If  IT9  =  1,  peaks  of  master  blade  flapping,  wing-tip  beamwise 
deflection,  wing-tip  chordwise  deflection,  and  first  normal 
coordinate  of  wing/pylon  are  selected  versus  time.  They  are 
printed  at  the  end  of  each  run. 

g .  Details  of  Rotor  and  Wing  Aerodynamics 

This  group  is  printed  at  the  end  of  the  second  rotor  revolu¬ 
tion  for  all  runs.  Instantaneous  aerodynamic  details  at  all 
blade  and  wing  stations  are  printed  for  reference. 

(1)  Rotor  Airloads  Details  (see  Figure  7) 

VP  is  the  vertical  component  of  the  blade 

local  velocity  (in/sec) 

VT  is  the  tangential  component  of  the 

blade  local  velocity  (in/sec) 


D  LIFT 


Figure  7.  Blade  Element  Aerodynamic 
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Details . 


PHI 

AL 

CL, CD 
CM 

MACH 

DHF 

D  Q 

DRADF 

DTHST 

DLIFT 

DDRAG 

DMOMT 

(  2 )  Wing  Airload 
VP 

VT 

PH 

AL 


is  the  ini  low  angle  (rad) 

is  the  angle  of  attack  (rad) 

is  the  local  blade  lift,  drag,  and 
pitching  moment  coefficients  respectively 

is  the  local  Mach  number 

is  the  rotor  in- plane  aerodynamic  drag 
force,  it  is  positive  if  it  tends  to  slow 
down  the  rotor  (lb) 

is  the  aerodynamic  pitching  moment  about 
the  blade  elastic  axis,  positive  for 
leading  edge  nose  up  (lb) 

is  the  blade  aerodynamic  radial  force  in 
shaft  plane.  Positive  if  it  points  to 
the  blade  root  (lb) 

is  the  blade  aerodynamic  thrust  perpen¬ 
dicular  to  shaft  plane.  Positive  if  it  is 
outward  from  the  shaft  plane  (lb) 

is  the  local  aerodynamic  lift  at  the 
quarter  chord.  It  is  perpendicular  to 
the  local  total  velocity  and  is  calculated 
directly  from  (lb) 

is  the  blade  local  aerodynamic  drag  at  the 
quarter  chord.  It  is  parallel  to  the 
local  total  velocity  and  is  calculated 
directly  from  Cq  (lb) 

is  the  blade  local  aerodynamic  pitching 
moment  at  the  quarter  chord.  It  is 
positive  if  blade  leading-edge  is  up,  and 
is  calculated  directly  from  (in- lb) 

Details  (see  Figure  8) 

is  the  vertical  component  of  the  wing 
local  velocity  (in/sec) 

is  the  tangential  component  of  the  wing 
local  velocity  (in/sec) 

is  the  local  inflow  angle  (rad) 

is  the  local  angle  of  attack  (rad) 
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Figure  8.  Wing  Element  Aerodynamic  Details. 


CL  is  the  local  lift  coefficient  for  the 

wing,  calculated  by 


where 

9L/3a  is  the  lift  curve  slope,  an  input 
on  Card  11  and,  a  is  the  local  angle  of 
attack 

CD  is  the  local  drag  coefficient 

WL  is  the  local  aerodynamic  lift  force 

WD  is  the  aerodynamic  drag  force 

WT  is  the  local  thrust  force  perpendicular 

to  the  free  stream  (lb) 

WH  is  the  wing  local  drag  force  parallel  to 

the  free  stream  (lb) 
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WQ  is  Lhe  local  torsional  moment  about  the 

elastic  axis,  positive  leading-edge  is  up 
( in- lb) 

The  machine  plots  are  sell-explanatory.  An  example  is  given 
in  Appendix  I . 

4  .  Suggested  User  Techniques 
a.  Data  Deck  Build-up 


In  preparing  a  data  deck  lor  program  ARAP06,  the  following 
steps  are  suggested. 

(1)  Start  With  Rotor  Rigid-Body  Freedom  -  The  user  should 
start  with  rotor  rigid-body  equations  TOPTION  IT5  =  0).  The 
data  required  for  this  option  consist  of  the  blades  inertia 
and  aerodynamic  properties,  i.e.  the  blade  twist,  mass  distri¬ 
bution,  mass  unbalance,  variation  of  blade  chord  with  blade 
radius,  rotor  preconc,  pitch-flap  coupling,  and  blade  collec¬ 
tive  pitch.  Tabulated  forms  for  the  blade  section  Cl,  Cq ,  and 
Cm  are  also  needed.  These  data  form  a  fundamental  data  deck 
for  the  program  which  can  be  readily  checked  out. 

(2)  Add  Wing  Elastic  Degrees-of-Freedom  -  After  checking  out 
the  rotor  rigid-body  degrees-of-freedom,  the  wing  flexible 
motion  should  be  added  (IT5  =  1).  This  requires  that  the  wing/ 
pylon  inertial  properties,  linear  dimensions  and  normal  modes 
as  well  as  natural  frequencies  be  supplied.  The  lift-curve 
slope  of  wing  airfoil  must  also  be  estimated.  The  wing 
equations  should  first  be  checked  out  with  the  proprotor  blade 
segment  lengths  set  to  zero  which  eliminate  the  proprotor 
aerodynamics.  Checkout  should  include  verification  of 
coupled  mode  shapes  and  frequencies. 

(3)  Add  Blade  Elastic  Degrees-of-Freedom  -  The  proprotor  blade 
elastic  motion  should  be  added  last.  This  requires  that  the 
blade  natural  frequencies,  normal  modes,  and  normalized  blade- 
bending  moments  for  each  mode  be  supplied  in  tabulated  forms 
(IT4  should  be  0). 

(4)  Other  Options  -  The  various  ARAP06  options  such  as  those 
for  wing/rotor  aerodynamic  interference,  blade  feathering  or 
gust  encounter  should  be  used  only  after  the  basic  math  model 
has  been  checked  out.  Many  of  the  options  require  hand 
calculation  prior  to  their  use.  For  example,  if  option  IT2 
is  to  be  used  to  feather  or  unfeather  the  proprotor,  the 
proper  blade  collective  pitches  should  be  calculated  by 
equation  134  in  Volume  I.  Likewise  for  the  1-cosine  gust  en¬ 
counter  option,  IT3  =  2,  the  gust  encounter  period  may  be  cal¬ 
culated  by 
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At 


25  WC 
U 


in  which  WC  is  the  wing  mean  aerodynamic  chord,  U  is  the 
velocity  along  flight  path.  These  are  input  parameters  on 
card  11  and  card  05,  respectively. 

b.  Machine  Time  Requirements 

The  ratio  of  the  machine  time  to  the  simulated  real  time  depends 
upon  flight  conditions  and  options.  Since  the  equations  of 
motion  in  ARAP06  are  nonlinear  in  nature  and  are  solved  simul¬ 
taneously  by  iterations,  the  more  complicated  the  input  con¬ 
dition  the  greater  the  number  of  iterations  required.  A  few 
examples  of  machine  time  requirements  are  given  below.  (The 
machine  time  here  referred  is  based  upon  an  IBM  360-65.) 

Example  1  -  The  blade  rigid-body  flapping  and  rotor  angular 
freedom  are  locked  out.  The  blades  are  inelastic  and  the 
proprotor  is  coupled  with  the  wing/pylon  degrees-of-freedom 
(4  modes).  Forty  seconds  of  machine  time  for  one  second  of 
real  time  was  required. 

Example  2  -  With  the  blade  flapping  and  rotor  angular  freedoms 
added,  the  machine  time  is  increased  to  ninety  seconds  per 
secound  of  real  time. 

Example  3  -  With  the  blade  rigid-body  flapping  and  rotor 
angular  freedoms  locked  out,  but  with  flexible  blades  (4  modes) 
coupled  with  wing/pylon  motions  (4  modes) ,  approximately  six 
minutes  of  machine  time  is  required  for  one  second  of  real  time. 

Example  4  -  If  the  proprotor  blade  flapping  and  angular  freedom 
is  added  to  example  3,  approximately  nine  minutes  of  machine 
time  is  required  for  one  second  of  real  time. 

Proper  choice  of  itegration  step  has  impact  on  both  the  machine 
time  and  program  accuracy.  Too  small  an  integration  step  in¬ 
creases  the  machine  time.  However,  too  large  integration  step 
will  prevent  the  program  from  giving  accurate  results.  The  user 
should  experiment  with  integration  step  and  establish  the 
largest,  accurate  step. 

c.  Acceleration  Averaging  Requirements 

Numerical  difficulties  are  occasionally  encountered  with  ARAP06 
because  of  the  nonlinearity  of  the  equations  of  motion.  This 
manifests  itself  as  a  higher  frequency  vibration  at  the  time 
step  frequency.  A  weighted  average  of  the  accelerations  can  be 
used  to  eliminate  this  problem.  This  option  is  controlled  by 
IT14. 
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B.  DFAL17  -  FPR  Blade  Folding  Flutter  Analysis 

The  inputs  to  DFAL17  consist  of  the  geometry,  mass  and  stiff¬ 
ness  distributions  of  the  wing,  pylon,  and  blades.  The  follow¬ 
ing  section  properties  can  vary  with  blade  fold  angle:  strip 
width,  chord,  reference  semispan  and  chord,  cosine  of  blade 
fold  angle,  mean  aerodynamic  chord,  distance  between  dumbbells, 
blade  sweep,  distance  from  elastic  axis  to  dumbbells,  mass 
unbalance,  and  inertia. 

The  card  order  in  a  DFAL17  data  deck  is  shown  in  Figure  9.  If 
the  program  is  being  used  to  calculate  natural  frequencies  and 
mode  shapes,  the  aerodynamic  data  set  is  not  required. 

1 .  Input  Format 

CARD  01  Format  20A4: 

Title  Card;  Column  1-80,  alphameric. 

CARD  02  Format  14: 

Col.  1  -  4  NSIZ  Matrix  size  =74,  if  minus, 

punch  flexibility  matrix 
and  eigenvalues 

Col.  5  -  8  NAERO  Reduced  velocities  (1/K) , 

if  NAERO  =  0,  skip  Cards 
3-19  (aero  input) 


Col. 

9  - 

12 

MODES 

Number  of  eigenvalues  to 
be  calculated 

Format  E12 

.5: 

Col. 

13  - 

24 

FK 

Flex  scaling  factor;  if 

FK  =  0  then  a  value  of 

FK  =  1  is  assumed 

Col. 

25  - 

36 

AMK 

Mass  scaling  factor; 
if  AMK  =0,  FM  =  1 

Col. 

37  - 

48 

RHO 

Density;  if  RHO  =  0,  (slug/ft3) 

RHO  =  0.002378 

Cards  03  through  19  contain  aerodynamic  inputs,  omit  Cards  03 


through 

19  if  NAERO  =  0. 

CARD  03 

Format  E12.5: 

Col. 

1-12 

BR 

Reference  semichord 

(ft) 

Col. 

13  -  24 

s 

Reference  semispan 

(ft) 

Col. 

25  -  36 

BREF 

Reference  semichord 

(ft) 

Col. 

37  -  48 

Reduced  velocities  in 

Col. 

49  -  60 

VELCTY 

fields  of  E12.5  until 

Col. 

61  -  72 

NAERO  values  have  been 

entered 
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Figure  9.  Card  Order  in  DFAL17  Data  Deck. 
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CARD  04 

CARD  05 

CARD  06 

CARD  07- 

CARD  08 

CARD  09 

CARD  10 

CARD  11 

Indicates 


Format  E12.5: 


Col.  1  -  12 

COSLMA 

Cosine  of  wing 
sweep  angle. 

Col.  13  -  24 

BRA 

Reference  wing  semichord 

(ft) 

Col.  25  -  36 

SA 

Reference  wing  semi  span 

(ft) 

Col.  37  -  48 

CAPA 

Reference  wing  area 

(  f  1 2  ) 

Col.  49  -  60 

Format  6E12.5 

CBAA 

Wing  mean  aerodynamic 
chord 

(ft) 

Col.  1-72 

Format  El 2. 5: 

Col.  1-12 

Format  6E12.5 

DELTA Yi 

Width  of  wing  panels 
(5  wing,  2  pylon) 

(ft) 

Col.  1-72 

Format  El 2. 5: 

Col.  1  -  12 

Format  6E12.5 

Bi 

• 

« 

Wing  panel  semichords 
(5  wing,  2  pylon) 

(ft) 

Col.  1-72 

Format  El 2. 5: 

Col.  1-12 

Di 

Distance  between  dumb¬ 
bells  on  wing  &  pylon 
(5  wing,  2  pylon) 

(ft) 

Format  E12.5: 

Col.  1-12 

cosLm 

Cosine  of  blade  fold 
angle 

Col.  13  -  24 

BRA 

Reference  semichord 

(ft) 

Col.  25  -  36 

SA 

Reference  semispan 

(ft) 

Col.  37  -  48 

CAPA 

Reference  area 

(ft2) 

Col.  49  -  60 

Format  5E12.5 

CBAA 

Mean  aerodynamic  chord 

(ft) 

Col.  1  -  60 

Format  5E12.5 

DELTA Yi 

• 

• 

Blade  1,  section  widths 

(ft) 

Col.  1  -  60 

Format  5E12.5 

Bi 

Blade  1,  section 
semichords 

(ft) 

Col.  1-60 

Di 

Blade  1,  distance 
between  dumbbells 

(ft) 

more  than  one  card  may  be  required  for  these  data. 


41 


CARD 

12 

Format  E12.5: 

Col.  1  -  12  COSLMA 

Col.  13  -  24  BRA 

Col.  25  -  36  SA 

Col.  37-48  CAPA 

Col.  49  -  60  CBAA 

Cosine  of  blade  fold 
angle 

Reference  semichord 
Reference  semispan 
Reference  area 

Mean  aerodynamic  chord 

(ft) 

(iV 

(ft2) 
( ft ) 

CARD 

13 

Format  5E12.5: 

Col.  1-60  DELTAY  ^ 

Blade  2,  section  widths 

(ft) 

CARD 

14 

Format  5E12.5: 

Col.  1-60  Bi 

Blade  2,  section 
semichords 

(ft) 

CARD 

15 

Format  5E12.5: 

Col.  1-60  Di 

Blade  2,  distance 
between  dumbbells 

(ft) 

CARD 

16 

Format  E12.5: 

Col.  1-12  COSLMA 

Col.  13  -  24  BRA 

Col.  25  -  36  SA 

Col.  37  -  48  CAPA 

Col.  49  -  60  CBAA 

Cosine  of  blade  sweep 
angle 

Reference  semichord 
Reference  semi span 
Reference  area 

Mean  aerodynamic  chord 

(ft) 

(ft) 

(ft2) 

(ft) 

CARD 

17 

Format  5E12.5: 

Col.  1-60  DELTAYi 

Blade  3,  section  widths 

(ft) 

CARD 

18 

Format  5E12.5: 

Col.  1-60  Bi 

Blade  3,  section 
semichords 

(ft) 

CARD 

19 

Format  5E12.5: 

Col.  1-60  Di 

Blade  3,  distance 
between  dumbbells 

(ft) 

Structural 

Input 

CARD 

20 

Format  El 2. 5: 

Col.  1-12  WGSWEP 

Col.  13  -  24  BLDSWP 

Wing  sweep,  forward  is 
minus 

Blade  sweep,  aft  is  plus 

(deg ) 

(deg) 
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CoL.  25  -  36 

VRPYLS 

Col.  37  -  48 

LTPYLS 

CoL.  49  -  60 

FOLDS  F 

Col.  61  -  72 

SMASTF 

CARD 

21 

Format  E12.5: 

Col.  1  -  12 

FXCLST 

Col.  13  -  24 

MVCLST 

Col.  25  -  36 

LCA 

Col.  37  -  48 
Col.  49  -  60 
Col.  61  -  72 

PYLOFF 

H 

EHUB 

CARD 

22 

Format  E12.5: 

Col.  1-12 

LW 

Col.  13  -  24 

EIBW 

Col.  25  -  36 

EICW 

Col.  37  -  48 

GJW 

CARD 

23 

Format  E12.5: 

Col.  1  -  12 

LW 

Col.  13  -  24 

EIBW 

Col.  25  -  36 

EICW 

Col.  37  -  48 

GJW 

CARD 

24 

Format  E12.5: 

Col.  1  -  12 

LW 

Col.  13  -  24 

EIBW 

Col.  25  -  36 

EICW 

Col.  37  -  48 

GJW 

Pylon  fore  and 
aft  stiffness 

( in- lb/ rad ) 

Pylon  lateral 
stiffness 

( in- lb/rad ) 

Blade  fold 

( in- lb/rad ) 

stiffness 

Mast  torsional 
stiffness 

( in- lb/rad) 

Fixed  control 

( in- lb/rad) 

stiffness 
Movable  control 
stiffness 

( in- lb/rad) 

Distance  elastic 
axis  to  conversion 

(  in) 

axis 

Pylon  offset 

(  in) 

Mast  length 

Hub  offset 

(in) 

Length  wing 

( in) 

segment  1 

El  beam  wing 

(lb-in2) 

segment  1 

El  chord  wing 

(lb-in2) 

segment  1 

CG  wing  segment 

1 

( lb- in2) 

Length  wing 

(in) 

segment  2 

El  beam  wing 

(lb- in2) 

segment  2 

El  chord  wing 

( lb- in2 ) 

segment  2 

GJ  wing  segment 

2 

( lb-in^ ) 

Length  wing 

(  in) 

segment  3 

El  beam  wing 

( lb-in2) 

segment  3 

o 

El  chord  wing 

(lb-in2) 

segment  3 

o  . 

GJ  wing  segment 

3 

( lb-m2 ) 
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CARD  25 


Format  EL2 

.5: 

CoL . 

1  - 

12 

LW 

Col . 

13  - 

24 

El  BW 

Col. 

25  - 

36 

EICW 

Col. 

37  - 

48 

G.JW 

CARD  26  Format  E12.5: 


Col. 

1  - 

12 

LW 

Col. 

13  - 

24 

El  BW 

Col. 

25  - 

36 

EICW 

Col. 

37  - 

48 

GJW 

Length  wing 
segment  4 

(  i  n ) 

( 1 b- in2  ) 

El  beam  wing 
segment  4 

El  chord  wing 
segment  4 

GJ  wing  segment  4 

(lb- in  ^ ) 

( Ib-i n2 ) 

Length  wing 

(in) 

segment  5 

El  beam  wing 

( lb- in2) 

segment  5 

El  chord  wing 

( lb-in2) 

segment  5 

GJ  wing  segment  5 

( lb-in2 ) 

CARD 

27 

Format  E12.5: 

Col.  1-12 

LB 

Col.  13  -  24 

EIBB 

Col.  25  -  36 

EIBC 

Col.  37  -  48 

GJB 

CARD 

28 

Format  E12.5: 

Col.  1-12 

LB 

Col.  13  -  24 

EIBB 

Col.  25  -  36 

EIEC 

Col.  37  -  48 

GJ 

CARD 

29 

Format  E12.5: 

Col.  1-12 

LB 

Col.  13  -  24 

EIBB 

Col.  25  -  36 

EIBC 

Col.  37  -  48 

GJB 

CARD 

30 

Format  El 2. 5: 

Col.  1  -  12 

LB 

Length  blade 

( in) 

segment  1 

El  beam  blade 

(lb- in2) 

segment  1 

El  chord  blade 

( lb-in2 ) 

segment  1 

GJ  blade  segment 

1 

( lb-in2 ) 

Length  blade 

( in) 

segment  2 

El  beam  blade 

(lb-in2) 

segment  2 

El  chord  blade 

( lb-in2 ) 

segment  2 

GJ  blade  segment 

2 

( lb-in2) 

Length  block 

(in) 

segment  3 

El  beam  blade 

( lb- in2) 

segment  3 

El  chord  blade 

(lb-in2) 

segment  3 

GJ  blade  segment 

3 

(lb-in2) 

Length  blade 

( in ) 

segment  4 
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CoL.  13  -  24  EIBB 

El  beam  blade 

( lb-ii.2) 

segment  4 

CoL.  25  -  36  EIBC 

El  c ho xxl  blade 
segment  4 

(lb-in**) 

rj 

CoL.  37  -  48  G.JB 

GJ  blade  segment  4 

(lb- in^ ) 

CARD  31 

Format  E12.5: 

CoL.  L  -  L2  LB 

Length  blade 
segment  5 

(  in ) 

CoL.  L 3  -  24  EIBB 

El  beam  blade 
segment  5 

(lb- in^ ) 

CoL.  25  -  36  EIBC 

El  chord  blade 
segment  5 

( lb-in^) 

CoL.  37  -  48  GJB 

GJ  blade  segment  5 

( lb-in^ ) 

Mass  Input 

CARDS 

Format  6EL2.5: 

32  &  33 

CoL.  L  -  72  X± 

Distance  EA  to  for¬ 
ward  central  point -- 
12  total,  5  wing,  2 
pylon,  5  blade  is 
positive 

(  in) 

CARDS 

Format  6EL2.5: 

34  &  35 

CoL.  L  -  72  Yi 

Distance  EA  to  aft 
central  point--12 
total,  5  wing,  2 

(in) 

pylon,  5  blade  is 
positive 

CARDS 

Format  6E12.5: 

36  &  37 

CoL.  L  -  72  M 

Weight  to  be  lumped 
at  each  dumbbell 

(lb) 

CARDS 

Format  6E12.5: 

38  &  39 

CoL.  L  -  72  S 

Unbalance  to  be 

( in- lb) 

lumped  at  dumbbell 
is  minus  if  forward 
of  EA 

CARDS 

Format  6E12.5: 

40  &  41 

Col.  1-72  I 

Inertia  to  be  lumped 
at  each  dumbbell 

( in-lb^ ) 
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2 .  Guide  to  Input  Format 


CARD  01 


CARD  02 


CARD  03* 


^Indicates 


TTL  is  used  to  identify  case  being  run. 

All  alphabetic  or  numeric  information  is 
reproduced  and  centered  on  first  output 
page. 

NSIZ  is  the  matrix  size  (74  max).  If  this 

number  is  negative  (-74),  then  the 
structural  influence  coefficient  matrix 
and  all  74  eigenvalues  arc  punched  on  cards 
(6E12.5  format)  for  use  as  input  to  DFAL18. 
The  flexibility  matrix  and  eigenvalues  are 
punched  for  the  first  1/K. 

NAERO  is  the  number  of  reduced  velocities  (1/K) 

<  15.  If  NAERO  =  0,  cards  3  through  19 
are  omitted.  If  NAERO  is  negative,  the 
aerodynamic  matrix  is  printed  out. 

MODES  is  the  number  of  eigenvalues  to  be  printed 
out  and  the  number  of  eigenvectors  to  be 
calculated. 

FK  is  the  flexibility  (SIC)  matrix  scaling 

factor.  The  El  and  GI  will  be  multiplied 
by  this  factor  to  reduce  input  data  require 
ments.  It  can  also  be  used  for  quick 
stiffness  changes  in  parametric  studies. 

If  FK  =  0  or  blank,  a  1  is  assumed. 

AMK  is  the  mass  matrix  scaling  factor.  The 

section  weights,  inertias,  and  unbalances 
are  multiplied  by  AMK.  This  can  be  used 
for  quick  mass  changes  in  parameter  changes 
If  AMK  =  0  or  blank,  a  1  is  assumed. 

RHO  is  the  air  density  (slug/ft^).  It  is  used 

in  the  calculation  of  airloads.  If  RHO  =  0 
then  0.002378  is  assumed. 

BR  is  the  reference  semichord  (ft).  The 

calculated  aerodynamics  are  multiplied  by 
BR  to  form  oscillatory  airloads. 

S  is  the  reference  semispan  (ft).  The 

calculated  aerodynamics  are  multiplied  by 
S  to  form  oscillatory  airloads. 

BREF  The  reference  semichord  (ft).  Usually 
the  same  as  BR.  Is  used  only  in  the 
velocity  calculation 

more  than  one  card  may  be  required  for  these  data. 
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BREF 


VELCTY 


CARD  04 


COSLMA 

BRA 

SA 

CAPA 

CBAA 


V  =  I 
v  K 


'R 


where  =  real  part  of  the  eigenvalue. 

V 

is  the  reduced  velocity  (1/K  =  -k — )  values. 

br" 

There  will  be  three  on  Card  3  and  six  on 
each  succeeding  card  until  NAERO  reduced 
velocities  arc  input. 

is  the  cosine  of  the  wing  sweep  angle. 

is  the  reference  wing  semichord  (ft). 

This  should  be  same  as  BR  on  Card  3. 


is  the  reference  wing  semispan  (ft) 
should  be  the  same  as  S  on  Card  03. 


This 


is  the  reference  area  (ft^).  This  should 
be  the  total  area  of  the  wing  and  three 
blades . 

is  the  wing  mean  aerodynamic  chord  (ft) . 

CBAA  =  |  cT  [  e  +  six] 

,  chord  root 

where  e  =  cHorTHp" 

CT  *  chord  tip 


CARD  05*  DELTAY^  is  delta  y  (the  strip  width  of  the  wing 

segments.  There  are  seven  strips--five 
for  the  wing  and  two  representing  the 
pylon.  Strip  six  is  for  pylon  pitch  air¬ 
loads  and  strip  seven  for  pylon  yawing 
airloads.  These  values  can  be  scaled  for 
3-D  lift  distribution  using  the  formula 


ai 

DY.  = 

1  2Ti 


UY. 


where  ai  is  the  estimated  lift  curve 
slope  of  the  individual  segment  (sweep 
effects  neglected),  DYi  is  the  actual 
geometric  strip  width. 

Two  cards  are  required  since  there  are 
seven  strips.  DYi  can  be  zero  if  desired. 
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CARD  06 


CARD  07 


CARDS  08, 
09,  10, 

LI 


CARDS  12, 
13,  14, 

15 

CARDS  16, 
17,  18, 

19 

CARD  20 


B. 

l 


D, 


( Blade 
1  aero 
data) 


( Blade 

2  aero 
data) 

(Blade 

3  aero 
data) 

WGSWEP 


BLDSWP 


VRPYLS 


LTPYLS 


FOLDSF 


SMASTF 


is  the  wing  section  seinichord  (it).  These 
are  the  semichords  for  each  wing  segment 
plus  the  two  pylon  sections.  Two  cards 
are  required.  must  be  nonzero. 

is  the  distance  (ft)  between  the  forward 
and  aft  dumbbell.  Normally  the  forward 
control  point  is  on  the  quarter  chordline 
and  the  aft  control  point  is  on  the  three- 
quarter  chordline.  This  makes  D^  equal  to 
B^.  D^  must  be  nonzero. 

same  format  as  cards  04  through  07  except 
blade  parameters  are  used.  There  are  only 
5  sections  for  each  blade,  so  there  are 
only  5  DY£,  5  ,  and  5  D^ . 

same  format  as  cards  08  through  11. 


same  format  as  cards  08  through  11. 


is  the  wing  sweep  angle  (degrees).  Forward 
sweep  is  negative. 

is  the  blade  fold  angle  (degrees).  Fold 
aft  is  positive.  Precone  is  a  negative 
blade  fold  angle. 

is  the  pylon  fore  and  aft  (pitch)  mounting 
spring  rate  (in-lb/rad).  This  number  is 
used  to  calculate  the  GJ  of  the  pylon 
offset  element. 

is  the  pylon  lateral  mounting  spring  rate 
(in- lb/rad).  This  stiffness  is  used  to 
calculate  the  chord  El  of  the  pylon  offset 
element . 

is  the  blade  fold  hinge  spring  rate 
(in-lb/rad).  This  stiffness  is  used  to 
calculate  the  EIC  of  the  one-inch  element 
located  between  the  fold  hinge  and  the 
first  blade  element. 

is  the  mast  torsional  stiffness  (in-lb/rad). 
This  number  is  used  to  calculate  the  mast 
GJ. 
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CARD  2  L 


FXCLST  is  the.  fixed  pitch  change  axis  spring  rate 
(in-lb/rad).  This  stiffness  is  used  to 
calculate  the  GJ  of  the  element  inboard 
of  the.  fold  hinge. 

MVCLST  is  the  moveable  pitch  change  axis  spring 
rate  (in-lb/rad).  This  stiffness  is  used 
to  calculate  the  GJ  of  the  one-inch 
element  located  between  the  blade  fold 
axis  and  the  first  blade  segment. 

LGA  is  the  distance  in  inches  between  the 

wing  elastic  axis  and  the  conversion  axis. 
Parallel  to  rotor  shaft.  Positive  if  con¬ 
version  axis  is  aft  of  the  elastic  axis. 

PYLOFF  is  the  lateral  distance  from  the  wing  tip 
to  the  shaft  centerline  (pylon  offset). 

H  is  the  distance  from  the  intersection  of 

the  conversion  axis  and  the  shaft  axis  to 
the  hub  (mast  length  -  in). 

EHUB  is  the  distance  from  the  intersection  of 

the  shaft  centerline  and  hub  to  the  blade 
fold  axis  (fold  hinge  offset  -  in). 

CARDS  contain  wing  structural  data.  Each  card  corresponds 

22-26  to  a  wing  segment.  The  format  for  each  card  is  the 

same. 

LW  is  the  length  of  ith  wing  segment  in  inches 

EIBW  is  the  beamwise  El  for  ith  wing  segment 
in  lb-in2 

EICW  is  the  chordwise  El  for  ith  wing  segment 
in  lb-in^ 

GJW  is  the  GJ  for  ith  wing  segment  in  lb-in2 

Each  wing  structural  segment  ends  at  a  lumped  mass 
point.  The  mass  point  should  be  located  in  the 
center  of  the  corresponding  aerodynamic  panel.  The 
distance  from  the  end  of  the  fifth  segment  to  the 
wing  tip  (segment  6)  is  assumed  to  be  one-half  the 
length  of  the  fifth  wing  segment.  The  stiffness 
of  segment  6  is  assumed  to  be  the  same  as  segment  5. 
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CARDS 

27-3L 


CARDS 

32-41 


contain  blade  structural  data.  Each  card  corresponds 
to  a  blade  segment.  The  blades  are  assumed  to  be 
identical.  The  format  ior  each  card  is  the  same. 

LB  is  the  length  of  ith  blade  segment  in 

inches 

EIBB  is  the  bcamwise  El  for  ith  blade  segment 

in  lb-mz 

EICB  is  the  chordwise  El  for  ith  blade  segment 

in  Ib-in^ 

GJB  is  the  GJ  for  ith  blade  segment  in  Ib-in^ 

Each  blade  segment  ends  at  a  mass  point  which  should 
be  located  at  the  center  of  the  corresponding  aero¬ 
dynamic  panel.  Therefore  the  fifth  blade  segment 
will  end  short  of  the  blade  tip.  The  blade  begins 
one-inch  outboard  of  the  fold  hinge. 

contain  inertial  data  for  the  wing  and  blades. 

X.  is  the  distance  in  inches  from  the  elastic 

axis  to  the  forward  control  point.  The 
forward  central  point  is  located  on  the 
quarter  chordline.  is  used  in  the 

generation  of  the  coupled  mass  matrix  and 
the  SIC  matrix.  X^  is  always  positive  and 
it  is  advisable  (for  matrix  conditioning) 
that  X-£  be  of  some  value  rather  than  zero. 
The  Xi  for  the  pylon  could  be  the  same  as 
the  mast  length  and  thus  show  the  hub 
motion.  The  five  wing  X^,  the  pylon  pitch, 
pylon  yaw,  and  five  X£  for  one  blade  are 
input  consecutively.  The  blades  are 
assumed  to  be  identical  so  the  blade  input 
is  not  repeated. 

Y.  is  the  distance  in  inches  from  the  elastic 

axis  to  the  aft  control  point.  The  aft 
control  point  is  normally  located  on  the 
three-quarter  chordline.  Yi  is  used  in 
the  generation  of  the  coupled  mass  matrix 
and  the  SIC  matrix.  is  used  in  the 

generation  of  the  coupled  mass  matrix  and 
the  SIC  matrix.  Y^  is  always  positive  and 
it  is  advisable  that  it  be  of  some  value 
rather  than  zero.  The  five  wing  Yi,  the 
pylon  pitch,  pylon  yaw,  and  the  Yi  for  one 
blade  are  input  consecutively.  The  blades 
are  assumed  to  be  identical  so  the  blade 
data  are  not  repeated. 
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M  is  the  section  weight  in  pounds.  This 

weight  is  distributed  on  the  fore  and  aft 
control  points  forming  a  dumbbell.  M 
cannot  be  zero.  M  for  the  five  wing 
segments,  pylon  pitch  and  yaw,  and  the 
five  blade  segments  are  in  consecutive 
order.  The  hub  weight  is  normally  lumped 
in  the  pylon  weight  rather  than  the  first 
blade  segment. 

S  is  the  section  unbalance  in  in-lb.  The 

unbalance  is  used  in  forming  the  mass 
matrix.  A  section  center  of  gravity  forward 
of  the  elastic  axis  produces  a  negative 
unbalance.  S  can  be  zero.  The  unbalance 
for  the  five  wing  segments,  pylon  pitch 
and  yaw,  and  the  five  blade  segments  are 
in  consecutive  order. 

I  is  the  section  inertia  I,  in  lb-in2.  This 

term  is  the  inertia  about  the  elastic  axis. 
This  parameter  cannot  be  zero.  The  inertias 
for  five  wing  segments,  pylon  pitch  inertia, 
pylon  yaw  inertia,  and  five  blade  segments 
are  entered  consecutively. 

3.  Output  Format  and  Guide 

The  output  of  DFAL17  consists  of  four  parts: 

-  All  input  data. 

-  Constants  which  are  internally  calculated.  These 
include  the  wing-tip  spring  rates,  the  flexibility 
and  mass  matrices,  and  the  aerodynamic  matrix,  if 
desired. 

-  Flutter  information  consisting  of  the  eigenvalues  and 
vectors  for  each  mode,  the  frequency  structural 
damping  required,  and  velocity  in  knots  associated 
with  each  frequency. 

-  Punched  cards  containing  the  structural  influence 
coefficient  matrix  and  eigenvalues,  if  desired. 

The  output  format  is  illustrated  by  a  sample  case  in 

Appendix  I. 
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vttNPp 


4 .  Suggested  User  Techniques 

a .  Problems  Encountered 

Any  large  computer  program  is  susceptible  to  numerical  in¬ 
accuracies.  DFAL17  is  no  exception.  Because  of  the  manipula¬ 
tion  of  Large  complex  matrices,  roundoff  errors  will  occur. 

The  numerical  difficulties  have  been  eliminated  as  much  as 
possible,  but  the  remaining  oddities  still  remain. 

The  tip  spring  rates  (beam  and  chord-- printed  output) 
will  be  different  for  various  blade  fold  angles.  This 
variation  is  small  (3-5  percent)  and  produces  little 
change  in  output  (frequency,  damping,  velocity,  and 
mode  shapes).  The  variation  is  caused  by  numerical 
inaccuracies  during  the  condensation  step. 

At  higher  reduced  velocities,  numerous  modes  will 
have  printed  "Imaginary  frequency"  in  lieu  of  a 
frequency.  Physically  this  signifies  a  divergence 
mode.  By  taking  smaller  reduced  velocity  increments 
prior  to  this  point,  the  velocities,  damping,  and 
frequencies  of  interest  can  be  calculated. 

All  data  must  be  checked  for  accuracy  since  input 
errors  may  produce  little  or  no  change  in  frequency 
but  gross  differences  in  damping  required  and  mode 
shapes . 

-  When  making  an  element  "rigid"  do  not  input  El  or  GJ 

more  than  10^  larger  then  the  stiffnesses  of  surrounding 
elements.  This  can  cause  extreme  numerical  errors. 

b.  Usage  Hints 

Given  below  are  some  hints  to  make  the  program  more  flexible. 

Aerodynamics,  as  calculated  by  the  program,  are 
2-D  and  incompressible.  By  variation  of  strip  widths, 
semichords,  and  density,  the  effect  of  3-D  theory, 
compressibility,  and  altitude  variations  can  be 
approximated.  (See  3-D  theory  example,  Figure  10.) 

Make  full  use  of  the  flexibility  and  mass  constants, 
the  spring  rate  inputs,  and  pylon  parameters.  These 
terms  will  enable  most  parameter  sweeps  to  be 
accomplished  by  single  number  change  each  run. 
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Example,  3-D  theory: 


The  surface  is  assumed  to  have  an  elliptic  lift  distribution 
as  shown  below: 


Figure  10.  3-D  Lift  Distribution. 


To  approximate  3-D  theory,  multiply  the  corresponding  AY's  by 
1.0,  0.95,  0.86,  0.73,  and  0.36. 

To  approximate  the  effect  of  compressibility,  multiply  the 

bi’s  by  the  ratio  of  %  compressible _  _ 

CL  incompressible 

The  altitude  can  be  approximated  by  varying  the  RHO  input  to 
correspond  to  the  altitude  desired. 
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c .  Time  Estimating 

Computer  time  for  this  program  is  a  function  of  modes  requested, 
time  sharing,  priority,  and  aerodynamics.  Generally,  for 
multiple  problems  within  a  run  allow  three  minutes  for  the 
first  case  and  two  minutes  for  each  additional  case.  Clock 
time  on  an  IBM  360/65  runs  about  2:1. 

C.  DFAL18  -  FPR  Blade  Folding  Response  Analysis 

Most  of  the  inputs  to  DFAL18  are  the  same  as  used  in  DFAL17. 

A  DFAL17  case  must  be  run  prior  to  running  DFAL18  to  generate 
the  flexibility  matrix  and  eigenvalues  inputs  to  DFALl 8 . 

The  card  order  in  a  DFALl 8  data  deck  is  shown  in  Figure  11. 


1.  Input  Format 
CARD  01  Format  2QA4: 

Title  Card;  Column 

CARD  02  Format  14: 

Col.  1  -  4  NG 


Col.  5  -  8  NT 


Col. 

9 

-  12 

NP 

Col. 

13 

-  16 

NSG 

1-80,  alphameric. 


If  NG  is  minus,  then 
all  response  points 
will  be  printed. 

If  NG  =  1,  sine  wave 
gust. 

If  NG  *  2,  sine- 
squared  gust. 

If  NG  *  3,  step  gust. 

If  NG  ■  4,  pulse  gust. 

If  NG  *  5,  step  gust 
with  Wagner  function. 

If  NT  *  0,  calculate 
coordinates  only. 

If  NT  =  1,  calculate 
force,  torque,  and 
moments  only. 

If  NT  =  2,  calculate 
coordinates  and  forces, 
torques,  and  moments. 

Number  of  input  flutter 
frequencies  NP  <  10. 

Number  of  coordinates 
to  be  run  NSG  <  10 . 
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Figure  11.  Card  Order  in  DFAL18  Data  Deck. 


Col.  17  -  20  NSCB 
21  -  27 


Coordinate  designa¬ 
tion  to  be  run, 

NSG  values. 


CARD  03  Format  E12.5: 


Col. 

1  - 

12 

RHO 

Col. 

13  - 

24 

V 

Col. 

25  - 

36 

CLA 

Col. 

37  - 

48 

VG 

Col. 

49  - 

60 

BR 

Col. 

61  - 

72 

S 

CARD  04  Format  E12.5: 


Col. 

1  - 

12 

DT 

Col. 

13  - 

24 

OM 

Col. 

25  - 

36 

TMAX 

Col. 

37  - 

48 

FK 

Col. 

49  - 

60 

WGSWEP 

Col. 

61  - 

72 

BLDSWP 

*  3 

Air  density  (slug/ft  ) 

Aircraft  forward  (knots) 

velocity 

Lift  curve  slope 

for  gust 

Gust  velocity  (ft/sec) 

Reference  semichord  (ft) 

Reference  semispan  (ft) 


Incremental  time  step  (sec) 

Blank,  gust  frequency  (cps) 

or  a  pulse  duration  (sec) 

Maximum  response  time  (sec) 

Flex  matrix  scale  factor 
Wing  sweep,  minus  (-)  (deg) 

forward 

Blade  sweep,  plus  (+)  (deg) 

aft 


Card  05  is  repeated  for  each  coordinate  to  be  analyzed,  i.e. 
NSG  cards  are  required. 


CARD  05 

Format  E12.5: 

Col.  1-12  EIB 

( lb-in^ ) 

Col.  13  -  24  GJ 

(lb-in^) 

Col.  25  -  36  Ll 

(in) 

Col.  37  -  48  L2 

(in) 

CARDS 

6-9 

Format  6E12.5:nv 

Col.  1  -  72 

(ft) 

CARDS 

10-13 

Format  6E12.5:R 

Col.  1  -  72 

(ft) 

CARD  14 

Format  E12.5: 

Col.  1-12  Gi 

Col.  13  -  24  coi 

(cps) 

Col.  25  -  36  G2 

Col.  37  -  48  Gl'2 

(cps) 

Col.  49  -  60  G3 

Col.  61  -  72  Cl 3 

(cps) 

Card  14  is  repeated  until  NP  damping  and  frequency  pairs  have 
been  input--three  to  a  card.  A  maximum  of  10  frequencies  is 
allowed . 
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Cards  15  through  30  are  identical  with  Cards  4  through  19 
of  DFAL17.  Cards  31  through  40  are  identical  with  Cards 
32  through  41  of  DFAL17. 


CARD  4L  Format  6E12.5: 

CoL.  1-72 

CARD  42  Format  6E12.5: 

Col.  1-72 

2 .  Guide  to  Input  Format 

CARD  01  TTL  is  the  title  card  consisting  of 

alphabetic  or  numeric  information 
required  for  identification  of  the 
run.  This  card  is  reproduced  on  the 
first  sheet  of  output  and  on  the 
plots . 

CARD  02  NG  is  a  control  number  designating  the 

gust  shape.  The  gust  begins  at 
t ime  =  0 . 


If  NG  >  0,  a  plot  tape  is  generated. 

If  NG  <0,  all  response  points  and 
loads  will  be  printed  as  well  as 
plotted. 

If  NG  =  1,  a  continuous  sine  wave 
with  frequency  defined  by  OM  (Card  4) 
in  cps  and  magnitude  defined  by  VG 
in  ft/sec. 

If  NG  =  2,  a  half  cycle,  1-cosine 
wave  with  frequency  defined  by  0M  in 
cps  and  magnitude  defined  by  VG 
in  ft/ sec. 

If  NG  =  3,  a  step  gust  input  with 
magnitude  defined  by  VG  in  ft/sec. 

If  NG  =4,  a  pulse  gust  with  a  dura¬ 
tion  defined  by  OM  in  seconds  with  a 
magnitude  defined  by  VG  in  ft/sec. 

If  NG  =  5,  a  step  gust  modified  by 
the  Wagner  function.  The  Wagner 
function  is  approximated  as: 
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L  -  0.L65  c-°*4555  _  0.335  e"0,3S 


where  S  = 

D 


V  =  forward  velocity,  ft/sec 

t  =  time  seconds 

B  =  reference  semichord,  ft 

The  magnitude  defined  by  VG  in  feet  per 
second . 


NT  is  an  option  designating  the  output 

required.  Consideration  should  be  given 
to  minimize  the  output  since  each  additional 
output  requires  computer  time. 

If  NT  =  0,  only  the  coordinate  responses 
will  be  calculated. 


NP 

NSG 

NSCP 


CARD  03  RHO 
V 

CLA 


VG 


If  NT  =  1,  only  the  forces,  moments,  and 
torques  will  be  calculated. 

If  NT  =  2,  the  forces,  moments,  torques, 
and  displacements  of  all  coordinates  will 
be  calculated. 

is  the  number  of  coupled  modes  to  be  input, 
NP  <  10. 

is  the  total  number  of  coordinates  to  be 
considered,  NSG  <  10. 

is  the  number  of  each  coordinate  to  be 
run.  There  will  be  NSG  coordinates 
entered  successively  in  14  format. 

is  the  air  density  in  slug/ft  . 

is  the  aircraft  forward  velocity  (knots). 

is  the  lift  curve  slope  of  the  lifting 
surfaces,  i.e.  wing  and  three  blades. 

This  lift  curve  slope  is  used  only  for 
the  gust  induced  airloads. 

is  the  maximum  vertical  gust  velocity  in 
feet  per  second. 
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CARD  04 


CARD  05 


BR 


S 


DT 


0M 


TMAX 


is  the  reference  semichord  (ft).  Should 
be  the  same  as  that  used  in  the  flutter 
analysis  DFAL17. 

is  the  reference  semispan  (ft).  Should  be 
the  same  as  that  used  in  the  flutter 
analysis  DFAL17. 

is  the  incremental  time  for  response  in 
seconds.  This  increment  should  be  small 
enough  to  allow  twelve  points  for  the 
smallest  frequency  of  interest.  Space  is 
provided  for  up  to  a  5-second  response  at 
0.0005  second  increments.  That  is,  1000 
time  points  can  be  run  for  each  coordinate. 
There  is  no  need  to  run  a  response  case 
for  long  periods  of  response  time  since 
the  points  of  interest  usually  are  the 
starting  transient  and  the  resultant  steady 
state.  Both  of  these  are  easily  achieved 
within  1000  time  points. 

depending  on  gust  option  (NG),  selected  0M 
could  be  a  blank  (step),  gust  frequency 
(sine  wave,  1-cosine,  cps)  or  a  pulse 
duration  (sec). 

Maximum  time  for  response  in  seconds. 


TMAX 

“DT“ 


<  1000 


FK  is  the  flexibility  matrix  scale  factor 

which  must  be  the  same  as  FK  used  in 
DFAL17  used  to  scale  the  flexibility  matrix. 

WGSWEP  is  the  wing  elastic  axis  sweep  angle  in 
degrees.  Forward  sweep  is  negative. 

BLDSWP  is  the  blade  elastic  axis  sweep  angle  in 
degrees.  Folding  aft  is  positive. 


One  card  must  be  supplied  for  each  coordinate  that 
is  to  be  analyzed.  This  card  provides  the  structural 
parameters  required  to  calculate  internal  loads  at 
the  coordinate. 


EIB  is  the  beamwise  El  for  the  segment 

corresponding  to  node  being  analyzed 
(Table  I)  in  lb-in2 
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GJ  is  the  torsional  stiffness  for  segment 

corresponding  to  node  being  analyzed, 

Table  I. 

L,  is  the  length  from  the  origin  to  the  left 

node  of  the  segment  being  analyzed.  Gan 
be  any  length  desired  but  must  be  /'L 
shorter  than  L2  ,  AL  =  L2  -  L]_.  For  the 
wing  L^,  is  the  length  from  the  root  to 
node  left.  For  the  blade  is  the  length 
from  the  intersection  of  the  mast  center- 
line  and  the  blade  elastic  axis. 

1^2  is  the  length  from  the  origin  to  the  right 

node  of  the  segment  being  analyzed.  Can 
be  any  length  desired,  but  must  be  AL 
longer  than  H,  AL  =  L2  -  Ll.  For  the  wing, 
Lg  is  the  length  from  the  root  to  the  node 
right.  For  the  blade,  L2  is  the  length 
from  the  intersection  of  the  mast  center- 
line  and  the  blade  elastic  axis. 

CARDS  DY  is  the  width  (in  feet)  of  the  aerodynamic 

6-9  surfaces  used  to  compute  the  gust  loadings 

on  the  wing  and  blades.  These  are  the  DYs 
used  in  DFAL17,  except  that  for  blades  2 
and  3,  the  DYs  must  be  one-half  of  those 
for  blade  1.  This  is  because  blades  2 
and  3  have  one-half  of  the  gust  velocity 
normal  to  their  planform.  These  DYs 
should  reflect  span  and  aspect  ratio 
effects.  The  DYs  are  entered  in  order  on 
four  succeeding  cards  in  the  order: 

CARD  06  Wing  -  five  segment  pylon  plot 

CARD  07  Pylon  yaw  and  blade  1  all  five  segments 

CARD  08  Blade  2  -  five  segments, 

Blade  3  -  one  segment 

CARD  09  Blade  3  -  segments  two  to  five. 

DY  can  be  zero. 

CARDS  B  is  the  local  semichord  (in  feet)  of  the 

10-13  aerodynamic  surfaces.  B  is  used  to  compute 

the  gust  loads  on  the  wing  and  blades. 

These  Bs  would  normally  correspond  to  the 
Bs  used  in  DFAL17.  The  Bs  are  entered  in 
order  on  four  succeeding  cards  in  the 
order: 
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CARD  10  Wing  (five  segments)  and  pylon  pitch 

CARD  11  Pylon  yaw  and  blade  1  (five  segments) 

CARD  12  Blade  2  (five  segments),  and  first  segment 
of  blade  3 


CARD  13  Blade  3  (segments  two  through  five). 

B  cannot  be  zero. 

CARD  14  NP  and  0)^  are  the  damping  factor  and 

frequency  in  cps,  respectively.  These 
are  obtained  from  the  DFAL17  output  for 
the  velocity  input  to  DFAL18.  NP  pairs 
must  be  input. 

Cards  15  through  30  are  identical  with  Cards  4  through  19  of 

DFAL17. 

Cards  31  through  40  are  identical  with  Cards  32  through  41 

of  DFAL17. 

CARD  41  FLEX  is  the  structural  influences  coefficient 

matrix  generated  with  DFAL17 .  This  matrix 
is  74  x  74  and  consists  of  962  cards, 
sequentially  numbered  and  in  6E12.5  format 
If  a  card  is  out  of  order,  the  program 
will  stop  and  print  an  error  message 
identifying  the  last  good  card. 

CARD  42  CMPLXR  are  the  complex  eigenvalues  for  the  first 

reduced  velocity  and  are  generated  using 
DFAL17.  The  reduced  frequency  should  be 
zero.  CMPLXR  consists  of  25  cards, 
sequentially  numbered. 

3.  Output  Format  and  Guide 


DFAL18  output  of  DFAL18  consists  of  four  parts: 

All  input  data. 

The  calculated  aerodynamic  matrix  and  eigenvalues. 

-  A  printed  time  history  of  the  coordinate  displacements 
forces,  moments,  and  torques. 

A  machine  (CALCOMP)  plot  of  displacements  and  the 
forces,  moments,  and  torques  at  the  right  and  left 
nodes . 


The  output  format  is  shown  by  a  sample  case  in  Appendix  I. 
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4. 


Suggested  User  Techniques 

a.  Problems  Encountered 

As  in  Program  DFALL7,  there  are  peculiarities  associated  with 
the  program.  Given  below  are  the  most  frequent  complaints. 

(1)  Loads  calculated  are  too  high  or  too  low.  This 
can  be  caused  by  any  of  the  following: 

(a)  Input  data  error;  card  out  of  place, 
integer  number  not  right  oriented  or 
floating  number  out  of  field. 

(b)  Incorrect  number  of  modes.  Since  blade 
modes  occur  in  threes,  to  include  one  or 
two  will  result  in  erroneous  loads. 

(c)  GJ  or  El  for  the  segment  is  wrong. 

(d)  Loads  are  correct  and  the  users  estimater  is  off. 

(2)  The  job  will  terminate  if  a  card  in  the  flexibility 
matrix  is  out  of  order.  All  cards  punched  as  output 
from  DFAL17  are  sequentially  numbered.  This  sequence 
number  is  checked  in  DFAL18. 

b.  Time  Estimate 

The  approximate  run  time  for  the  loads  at  one  node  is  approxi¬ 
mately  six  minutes.  This  time  will  vary  with  priority  time 
sharing  procedures  and  whether  output  is  printed  on  line  or 
off  line. 


SECTION  III 


PROGRAMMING  GUIDES 


This  section  contains  the  fundamental  structure  of  programs 
developed  under  this  contract.  Details  of  numerical  methods 
and  program  organization  are  then  discussed.  FORTRAN  listings 
will  not  be  included  in  this  section;  they  are  listed  in 
Appendix  II. 

A.  ARAP06 

I .  Method  of  Digital  Computer  Solution 

A  predictor-corrector  integration  technique  is  used  in  the 
solution  of  the  simultaneous  equations  of  motion  which  are 
nonlinear,  open-form  in  nature.  In  general  the  solution  of  a 
set  of  differential  equations  by  this  method  requires  the 
isolation  of  the  highest  order  derivatives.  For  example 
given  the  equation 


ax  +  bx  +  cx  =  F(t) 

where  the  coefficients  a,  b,  c  and  forcing  function,  F(t), 
may  all  be  time  dependent  and/or  nonlinear. 

To  solve  for  x,  the  equation  is  rearranged  in  the  following 
form 


Prior  to  evaluating  this  equation,  a  value  of  x,  called  iig ,  is 
guessed.  When  an  iterative  procedure  starts  at  a  new  time 
point,  Rg  is  used  from  the  solution  at  the  previous  time  point. 
Otherwise  it  is  predicted  by  use  of  point- slope  computational 
methods  to  minimize  the  difference  between  the  value  calculated 
and  the  value  guessed. 

Since  values  of  velocities  and  displacements  are  also 
required,  to  evaluate  the  equation  Taylor’s  series  is  used 
to  provide  a  numerical  integration  form  for  determining 
them  as  functions  of  the  guessed  acceleration. 

By  definition  the  integral  of  a  function  f'(x)  over  the 
interval  h  is : 


f (x+h) 


f  ( x) 


f ' (x)dx 
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Also  the  incremental  form  of  Taylor's  series  expansion  oi  a 
fraction  f(x)  is: 

2  3 

f  (x  +  h)  =  f  ( x )  +  hf'(x)  +  jj-f'Cx)  +  -^-f"'(x)  +  ... 

Thus  the  integral 

*+h  h2  k3 

f  '  (x)dx  =  hf  '  (x)  +  -jj-f  11  (x)  +  -^-f  '  '  '  (x)  +  ... 

x 

Defining  derivatives  of  f(x)  in  the  following  way: 

f ' (x)  =  f(x) 

We  can  produce  an  integration  form  which  permits  evaluation  of 
f(x+h)  in  terms  of  its  derivatives,  i.e.  numerical  integration. 
The  derivative  can  be  defined  in  other  ways,  for  example 

x:,  _  f (x+h)  -  f (x-h) 

f  (x) - ^ - 

But  for  simplicity  of  programming  the  first  definition  is 
used  here. 

Substituting  and  letting  h  =  At: 

f(x+At)  =  f(x)  +  At*f'(x)  +  — [f  '  (x+At )  -  f'(x)j 
2 

+  [f  .  »  (X+At)  -  f »  » (X)j  +  .  .  . 


Thus 


f(x+At)  =  f(x)  +  ^  f'(x+At)  +  f'(x)l 


Let 


+  [f  ’  '(x+At)  -  V  1  (  x )]  +  .  .  . 


f(x+At)  =  x 

O 
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retain  the  second  term  only,  one  obtains, 


Then  let 


x 


g 


X  ,  + 

n- 1 


+  X 

n- 


f(x+At)  =  x 

s 


one  obtains 


x 


g 


X 


n- 


I 


+ 


2 

\(/t)Z 

i)~r~ 


where  the  subscript  n-L  refers  to  the  previous  time  point. 

A  computed  acceleration  is  then  obtained  by  using  the  guessed 

values  of  x,  and  x  ,  i.c. 
g  g 


x  =  )  b  ^  £  x 

xc  a  a  xg  ‘  a  xg 

This  computed  value  in  general  will  not  be  equal  to  the 
guessed  one.  A  convergence  criterion  has  to  be  met  such  that 


where  t  is  the  error  tolerance. 

If  the  error  percentage  fails  to  be  less  than  e,  a  quadratic 
point-slope  subroutine  is  called  using  parabolic  curve  fitting 
to  minimize  the  error  percentage  to  zero.  A  better  predicted 
xg  is  then  obtained  for  the  next  iteration.  The  process 
repeats  itself  until  the  convergence  criterion  is  met. 
Throughout  these  prediction  and  correction  procedures,  the 
independent  variable-time  remains  constant.  When  the  error 
percentage  is  within  the  tolerance,  the  convergence  is  com¬ 
pleted.  The  value  of  time  is  then  increased  by  an  increment 
and  the  iterative  process  is  repeated. 

2 .  Subroutine  Descriptions 

a .  FORTRAN  Subroutines 

The  FORTRAN  subroutines  contained  in  the  Proprotor  Aeroelastic 
Analysis,  ARAP06 ,  are  listed  in  alphabetical  order  except  for 
the  main  program  which  is  first.  A  few  remarks  are  made  for 
each  subroutine,  which  indicate  the  general  purpose  or  use  of 
tne  subroutine  in  the  general  structure  of  the  program.  Great- 
detail  of  explanation  is  not  attempted  here  since  it  would  be 
redundant . 
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(1)  MAIN  -  This  routine  is  the  main  program  of  ARAP06.  This 
subprogram  sets  up  the  time  loop  and  directs  the  flow  of 
beginning  and  ending  of  the  whole  program.  The  rate  of 
feathering  or  unfeathering  of  rotor-blades  is  computed  here. 

(2)  AIRLDS  -  This  subroutine  handles  the  computation  of  rotor 
aerodynamic s . 

(3)  BEXIOM  -  Tabulated  forms  of  natural  frequencies,  normal 
modes  and  bending  moments  of  flexible  blades  as  functions  of 
rotor  rpm  and  blade  feathering  are  read  in.  Printing  out  of 
these  tables  is  also  included  in  this  subroutine. 

(4)  BLOCK  DATA  -  This  subroutine  contains  tabulated  forms 

of  blade  aerodynamic  lift,  drag,  and  moment  coefficients  Cl, 

Cp ,  and  Cm  as  functions  of  Mach  number  and  angle-of-attack. 

It  also  contains  the  tabulated  coefficients  of  blade  pitch¬ 
damping  as  functions  of  angle-of-attack  and  reduced  velocity. 

(5)  CLCDR  -  A  maximum  of  four  different  blade  airfoil  sections 
can  be  handled  in  ARAP06.  This  subroutine  determines  the  air¬ 
foil  for  each  blade  segment.  This  subroutine  also  provides 
the  interpolation  for  values  of  Cl,  Cp,  and  for  the  most 
inboard  blade  airfoil  section. 

li 

(6)  FLAPN1  -  Forcing  functions  for  blade  rigid-body  flappings 
and  rotor  rotational  degrees-of-freedom  are  calculated.  These 
equations  are  then  solved  simultaneously. 

(7)  FLAPN3  -  Tests  of  converging  criteria  for  blade  rigid-body 
flappings  and  rotor  rpm  are  conducted.  If  convergence  is 
completed,  forces  and  moments  due  to  rotor  aerodynamics,  rotor 
inertials  and  rotor  gravity  are  summed  and  ready  to  be  used  as 
a  forcing  function  for  wing/pylon  dynamic  equations.  Print¬ 
out  of  converged  rotor  variables  is  then  processed.  The  peaks 
of  master  blade  flapping  are  stored  if  required. 

(8)  FLEXBD  -  This  subroutine  solves  for  generalized  coordinates 
for  the  flexible  blades. 

(9)  FLEX1  -  Because  of  limitations  on  the  size  of  subprograms 
in  the  FORTRAN  compiler,  all  dynamic  terms  used  in  the  dynamic 
equations  of  motion  for  rotor-blades  cannot  be  handled  in  any 
single  subroutine.  Some  terms  are,  therefore,  calculated 

in  this  subroutine  and  will  be  used  in  both  subprograms 
FLAPN1  and  FLEXBD. 

(10)  FORGES  -  Rotor  inertial  loads  and  moments  are  calculated 
when  the  convergence  of  blade  generalized  coordinates  is 
completed. 
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(11)  HEADNG  -  This  subroutine  handles  the  FORMAT  of  the  print¬ 
out  headings .  The  input  options  have  control  on  the  selection 
of  variables  to  be  computed  as  well  as  on  the  output  headings. 

(12)  INPUT  -  The  read-and-write  of  input  data  is  done  in  sub¬ 
routines  BEXIOM  and  READWR.  The  subroutine  INPUT  changes 
input  parameter  units  if  necessary,  and  computes  time- 
invariant  quantities  from  input  data. 

(13)  FLTS  -  This  subroutine  plots  the  specified  variables. 

(14)  PSLOP  -  This  subroutine  minimizes  the  difference  between 
the  calculated  acceleration  and  the  predicted  one  and  provides 
a  better  prediction  for  next  iteration.  Let  x'g  be  the  pre¬ 
dicted  acceleration,  xc  be  what  calculated,  ana  Ax  =  xc  -  x'g 

a  point,  ( xg ,  Ax)  is  then  obtained  in  a  two-dimensional  space. 
When  iteration  loop  starts  at  a  new  time  point,  there  is  only 
one  (Rg,  AX)  available.  The  next  prediction  is  given  by 
simply  averaging  and  *XC  .  After  the  second  iteration,  there 
are  two  points  of  °(ttg,  Ax)  available.  A  straight  line  passing 
through  these  two  points  and  intersecting  the  axis  of  abscissa 
gives  a  value  of  Xg.  This  Kg  is  used  as  a  predicted  value  for 
the  next  iteration.  After  tne  third  iteration,  there  should 
always  be  three  points  of  (jig,  AX)  available,  a  quadratic 
equation  is  then  used.  The  Intersection  of  the  curve  repre¬ 
sented  by  this  quadratic  equation  and  the  abscissa  provides 
the  next  predicted  Xg  for  use  in  the  following  iteration.  Two 
subroutines  are  generally  obtained  by  solving  a  quadratic 
equation,  the  one  which  is  closest  to  the  last  guess  is  used. 

(15)  READWR  -  This  subroutine  reads  and  writes  all  the  input 
data  for  proprotor  and  wing/pylon  except  those  for  the  rctor 
flexible  blades. 

(16)  RSTN  -  Flapping  restraint  and  stops  of  a  gimballed 
rotor  are  handled  in  this  subroutine. 

(17)  TABINA  -  Interpolation  of  tabulated  forms  of  natural 
frequencies ,  normal  modes  and  bending  moments  of  flexible 
blades  as  functions  of  instantaneous  rotor  rpm  and  blade 
feathering  is  done  in  this  subroutine. 

(18)  TABINT  -  Excluding  the  most  inboard  blade  airfoil  section, 
the  Cl,  Cq,  and  Cm  of  the  rest  of  the  three  sections  of  a  blade 
are  provided  by  interpolation. 

(19)  WING  -  Wing  aerodynamics  and  generalized  force  for  wing/ 
pylon  system  are  calculated.  Dynamic  equations  of  motion  of 
wing/pylon  are  then  solved.  Converging  criteria  for  the 
solved  normal  coordinates  are  tested.  Pylon  accelerations 
after  completion  of  the  convergence  are  computed. 
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b.  Subroutines  Cross -Reference 

The  information  provided  in  Table  I  is  vital  lor  understanding 
of  the  subroutine  organization  of  program  ARAP06.  The  first 
column  contains  the  names  of  the  subroutines  in  alphabetical 
order,  except  for  the  main  program  which  is  first.  The 
second  column,  length,  contains  the  size,  in  bytes  in  hexa¬ 
decimal  number  base,  of  the  compiled  subprogram.  The  FORTRAN 
subroutines,  all  but  BLOCK  DATA,  have  been  compiled  on  the 
IBM  operation  system  360  FORTRAN  IV  (H)  Compiler  with  OPT  =  2. 
The  subroutine,  BLOCK  DATA,  has  to  be  compiled  on  FORTRAN  IV 
(G)  Compiler  because  of  its  size.  The  third  column  contains 
two  operations,  and  the  fourth  column  contains  the  names  of 
the  subroutines  associated  with  that  operation  and  given  in 
the  order  of  call  sequence. 

3 .  Flow  Outline 

The  flow  charts  in  Figures  13  and  14  describe  the  functional 
structure  of  the  program  and  the  logic  used  in  the  subroutine 
AIRLDS  respectively,  without  regard  to  flow  by  subroutines. 
Figure  13  provides  a  clear  picture  of  the  numerical  technique 
used  in  ARAF06  to  solve  the  equations-of-motion  simultaneously. 
Figure  14  demonstrates  how  the  wing-rotor  aerodynamic  inter¬ 
ference  is  handled  in  the  calculation  of  blade  airloads. 

Tedious  flow  charts  are  not  given  in  this  report.  The  reader 
can  refer  to  FORTRAN  listings  for  details. 
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TABLE  I 

ARAPOfa  SUBROUTINE  CROSS-REFERENCE 


Name 

Length 

Operation 

C r o  s s - Re ie  re  nc e 

MAIN 

CEO 

Calls 

INPUT,  HEADNG,  FLEXI, 
AIRLDS ,  FLEXBD  ,  FLAPN1  , 
FLAPN3  ,  PLTS 

AIRLDS 

147A 

Calls 

Is  Called 

By 

CLCDR 

MAIN 

BEXIOM 

1274 

Is  Called 

By 

READWR 

BLOCK 

DATA 

2A08 

Is  Used  By 

TABINT 

CLCDR 

930 

Calls 

Is  Called 

By 

TAB  I  NT 

AIRLDS 

FLAPN1 

1274 

Calls 

Is  Called 

By 

RSTN 

MAIN 

FLAPN3 

1814 

Calls 

Is  Called 

By 

PSLOP,  TABINA,  FORCES, 
RSTN,  WING,  HEADNG 

MAIN 

FLEXBD 

1ECA 

Calls 

Is  Called 

By 

PSLOP 

MAIN 

FLEXI 

94C 

.Is  Called 

By 

MAIN 

FORCES 

FIE 

Is  Called 

By 

FLAPN3 

HEADNG 

56C 

Is  Called 

By 

MAIN,  FLAPN3 

INPUT 

13A6 

Calls 

Is  Called 

By 

READWR,  TABINA 

MAIN 

PLTS 

2A2C 

Is  Called 

By 

MAIN 

PSLOP 

4  LA 

Is  Called 

By 

FLAFN3,  FLEXBD,  WING 

READWR 

24E6 

Calls 

Is  Called 

By 

BEXIOM 

INPUT 

RSTN 

304 

Is  Called 

By 

FLAPN1 ,  FLAPN3 

TABINA 

5F0 

Is  Called 

By 

INPUT,  FLAPN3 

TAB  I  NT 

9C2 

Is  Called 

By 

CLCDR 

WING 

12C0 

Calls 

Is  Called 

By 

PSLOP 

FLAPN3 
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Figure  12.  Flow  Chart  for  Program  Structure  -  ARAP06 . 
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Figure  12.  Concluded 
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Figure  13.  Flow  Chart  for  Subroutine  AIRLDS. 
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B. 


DFAL17 


1 .  Method  of  Solution 

Computer  program  DFAL17  solves  for  the  flutter  speed  of  a 
folding  proprotor.  The  stiffness  matrix,  mass  matrix,  and 
aerodynamic  matrices  are  constructed,  arranged  and  solved  with 
a  minimum  of  input  data.  The  user  needs  only  to  specify 
stiffnesses,  geometric  sizings,  and  fold  angles. 

The  dynamics  matrix  is  solved  by  ALLMAT,  a  QR  algorithum, 
which  solves  for  all  74  eigenvalues  in  a  minimum  of  computer 
time.  The  eigenvalues  are  then  arranged  in  order  of  decreas¬ 
ing  real  part  (increasing  frequency)  and  the  appropriate 
number  of  eigenvectors  are  calculated.  The  structural  in¬ 
fluence  coefficient  matrix  and  74  eigenvalues  can  be  punched 
output  for  use  as  input  to  computer  program  DFAL18. 

The  flow  of  the  computer  program  is  shown  by  a  macroscopic 
flow  diagram  (Figure  13).  A  brief  description  of  each  sub¬ 
routine,  a  flow  outline,  and  a  list  of  variable  arrays  are 
given  on  the  following  pages.  The  main  variable  arrays  are 
defined  at  the  end  of  this  chapter.  The  reader  is  referred 
to  the  FORTRAN  listing  (Volume  III)  for  program  details. 

2 .  Subroutine  Description 


SUBROUTINE  PRT3 

Initializes  and  builds  the  dynamics  matrix  using  an  assumed 
reduced  velocity  or  flutter  frequency. 

SUBROUTINE  PRT4 

Solves  the  dynamic  system  for  all  eigenvalues  and  requested 
eigenvectors.  Orders  eigenvalues  for  lowest  frequency,  out¬ 
puts  eigenvalues,  eigenvectors,  natural  frequencies,  and 
punched  eigenvalues  if  necessary. 

SUBROUTINE  DFAE ,  BESJ,  BESY 

This  subroutine  and  its  associated  routines  BESJ,  BESY,  form 
the  aerodynamic  influence  coefficient  matrix  corresponding  to 
a  given  reduced  velocity  or  flutter  frequency.  BESJ  computes 
J  Bessel  functions.  Methods  used  are  those  given  in  IBM 
SYSTEM/360  SCIENTIFIC  SUBROUTINE  PACKAGE,  VERSION  III 
PROGRAMMER’S  MANUAL,  H20-0205-3,  page  363-364. 


Fiyure  14.  DFAL17  Flow  Diagram. 
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SUBROUTINE  STGN 


Reads  stiffness  inputs,  forms  stiffness  matrix  parameters, 
inverts,  transforms,  and  produces  flexibility  matrix.  Punches 
flex  matrix  if  necessary. 

SUBROUTINE  ROTATE 

Builds  Large  stiffness  matrix,  180th  order,  and  transformation 
matrix. 

SUBROUTINE  MPRINT 

Prints  a  program  matrix  in  column  form. 

SUBROUTINE  REDUCE 

Eliminates  undesired  coordinates  from  the  symmetric  stiffness 
matrix  by  Gaussian  pivotal  reduction,  reduction  of  degrees  of 
freedom  from  180  to  74. 

SUBROUTINE  MULT 

Forms  a  matrix  product  UA  where  A  is  a  real  matrix,  in  this 
program  a  flex  matrix,  and  U  is  a  complex  aerodynamic-inertial 
matrix. 


SUBROUTINE  ALLMAT 

This  program  uses  a  QR  algorithm  with  special  attention  to 
minimization  of  required  storage  to  solve  for  the  eigenvalues 
and  eigenvectors  of  a  general  system  of  the  form  A  -  >1  =  0, 
where  A  is  a  complex  valued  matrix.  In  this  case  the  matrix 
A  is  74th  order,  and  not  symmetric. 


SUBROUTINE  INVERS 

This  routine  inverts  an  N  x  N  matrix  by  Gaussian  row  and 
column  pivot  operations. 


SUBROUTINE  TPROD 


T 

This  subroutine  performs  a  specific  matrix  triple ^ product  X  AX 
where  A  is  a  180th  order  symmetric  stiffness  matrix  in  upper 
triangular  form  and  X  is  a  quinary  diagonal  transformation 
matrix. 


SUBROUTINE  DPROD 


T 

terforms  a  matrix  triple  product  X  AX  where  A  is  a  symmetric 
N  x  N  matrix  and  X  is  a  general  transformation  matrix. 
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SUBROUTINE  ERRSET 


This  IBM  360  System  subroutine  suppresses  underflow  error 
print  -ats  which  occur  from  solution  of  an  undamped  dynamics 
matric  problem  with  an  eigenvalue  routine,  ALLMAT  which  is 
intended  for  solution  of  damped  systems.  In  some  cases,  the 
imaginary  part  of  the  eigenvalues  approaches  zero  slowly. 

3.  Flow  Outline 

MAIN 

(1)  Begin  reading  input  for  new  case.  Define  default 
parameters.  If  no  aero  input  is  to  be  called  for, 
to  to  (3). 

(2)  Read  aero  input  and  call  subroutine  DFAE  to  set  up 
aero  coefficients. 

(3)  Call  subroutine  STGN  to  compute  flexibility  matrix. 

(4)  Output  flexibility  matrix  and  scale  by  k  factor. 

(5)  Start  loop  to  evaluate  distributed  mass  matrix. 

(6)  Form  mass  matrix  by  triple  product  coordinate  trans¬ 
formation.  When  all  mass  elements  are  complete,  go  to 

(7). 

(7)  Print  mass  matrix  and  call  subroutine  PRT3  for  eigenvalue- 
eigenvector  computations. 

SUBROUTINE  PRT3 

(1)  Set  up  loop  to  build  and  evaluate  system  for  each 
reduced  velocity  input  in  MAIN  (2)  or  once  for  the 
structural  system  if  no  aero  exists. 

(2)  Zero  and  build  aero  matrix  transpose  by  calling 
subroutine  DF67. 

(3)  Add  the  mass  matrix  transpose  to  the  aero  matrix.  Units 
of  slugs  are  consistent. 

(4)  Form  product  (M+A)K-i,  the  augmented  mass  matrix  post- 
multiplied  by  the  flexibility  matrix. 

(5)  Transpose  the  dynamics  matrix  transpose  to  form  the 
dynamics  matrix. 

(6)  Call  subroutine  PRT4  to  compute  eigenvalues  and  eigen¬ 
vectors  of  dynamics  matrix. 
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(7)  If  aLL  reduced  velocities  have  been  run  return  to  MAIN, 
otherwise  go  back  to  (2)  for  next  velocity. 

SUBROUTINE  PRT4 

(1)  Compute  all  74  eigenvalues  by  calling  subroutine  ALLMAT. 

(2)  Use  loop  to  order  the  first  MODES  eigenvalues  in  ascending 
order  of  magnitude  of  the  real  or  frequency  component 
only.  This  corresponds  to  lowest  flutter  frequencies 
first . 

(3)  Print  eigenvalues.  If  punch  option  is  on,  punch 
eigenvalues . 

(4)  Start  loop  to  evaluate  eigenvectors. 

(5)  Call  subroutine  entry  ALLVEC  to  compute  the  eigenvector 
corresponding  to  a  given  eigenvalue. 

(6)  Normalize  the  eigenvector  by  division  by  the  largest 
element . 

(7)  Go  through  logic  to  print  eigenvectors  in  groups  of 
threes . 

(8)  If  all  eigenvectors  needed  have  been  computed,  go  to 

(9),  otherwise  to  (4)  for  next  eigenvector. 

(9)  Start  loop  to  e.valuate  all  frequency  and  damping  factors 
derived  from  flutter  frequencies  of  interest. 

(10)  Compute  and  print  flutter  frequency,  damping  and 
equivalent  velocity. 

(11)  At  end  of  loop  (9)  return  to  subroutine  PRT3, 

SUBROUTINE  DFAE 

(1)  Read  input  aero  data  case  and  print  all  quantities  input. 
Return  to  MAIN  program. 

(2)  At  entry  DF67  the  reduced  velocity  EKR  to  be  used  is 
input  to  the  subprogram. 

(3)  Start  loop  over  all  surfaces,  in  this  program,  four. 

(4)  Start  inner  loop  over  all  strips  on  each  surface. 

(5)  Set  up  transformation  matrix  A  dependent  on  surface-strip 
geometry. 
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(6)  Compute  J  and  Y  Bessel  functions  for  the  modified 
frequency . 

(7)  Set  up  the  oscillatory  aerodynamic  coefficient  matrix  B 
using  Bessel  functions  from  (6). 

T 

(8)  Do  a  transformation  matrix  triple  product  A  BA  to  get 
aerodynamic  influence  coefficient  matrix. 

(9)  Scale  the  coefficients  by  an  eiquivalent  air  mass 
factor. 

(10)  Use  a  loop  to  set  up  the  transpose  of  the  aerodynamic 
matrix  by  distributing  the  aerodynamic  influence 
coefficients  along  the  diagonal. 

(11)  Have  all  strips  been  computed?  If  not  return  to  (4). 

(12)  Have  all  surfaces  been  computed?  If  not  return  to  (3). 

(13)  Return  to  subroutine  PRT3  with  complex  aerodynamic 
matrix . 

SUBROUTINE  STG N 

(1)  Read  structural  inputs  and  print. 

(2)  Use  loop  to  initialize  angles  and  lengths. 

(3)  Set  up  quantities  required  for  stiffness  matrix. 

(4)  Gall  subroutine  ROTATE  to  build  180th  order  stiffness 
matrix  A  and  transformation  matrix  B. 

T 

(5)  Gall  subroutine  TPROD  to  form  triple  product  B  AB. 

(6)  Gall  subroutine  REDUCE  to  reduce  large  stiffness  matrix 
to  74th  order  by  elimination  of  coordinates. 

(7)  Read  parameters  required  for  mass  matrix  and  coordinate 
t  rans  f  orma  t io  n . 

(8)  Interchange  columns  11,  12,  and  13  in  stiffness  matrix 
and  initialize  transformation  matrix. 

(9)  Invert  stiffness  matrix  to  get  74th  order  flexibility 
matrix. 

(10)  Print  computed  stiffnesses  and  input  from  (7). 

(11)  Build  transformation  matrix  from  input  of  (7). 
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(12)  Call  subroutine  DPROD  to  do  transformation  triple 
product . 

(13)  If  punch  option  is  on,  punch  the  final  flexibility 
matrix . 

(14)  Return  to  MAIN. 


SUBROUTINE  ROTATE 

(1)  Zero  upper  triangular  stiffness  matrix  storage. 

(2)  Start  loop  to  build  180th  order  upper  triangular 
stiffness  matrix. 

(3)  Look  at  each  node  point,  30  in  all;  build  the  elements 
of  the  stiffness  matrix  by  coordinate  transformations  of 
a  coupled  2  node  stiffness  system. 

(4)  Store  all  computed  values  of  node  stiffness  quantities 
in  an  upper  triangular  form  due  to  symmetry. 

(5)  If  all  nodes  have  been  completed  go  to  (6),  otherwise 
to  ( 2)  . 

(6)  Initialize  node  independent  values  of  the  gamma  trans¬ 
formation  matrix. 

(7)  Build  the  gamma  transformation  matrix  as  a  set  of  3  x  3 
matrices,  2  for  each  node. 

(8)  When  the  entire  gamma  transformation  matrix  is  complete, 
return  to  subroutine  STGN. 

SUBROUTINE  REDUCE 

(1)  Start  loop  to  evaluate  all  180  degrees  of  freedom. 

(2)  If  the  selected  degree  of  freedom  is  not  to  be  reduced 

to  (8). 

(3)  Select  pivot  element  of  coordinate  to  be  reduced. 

(4)  Use  loop  to  build  a  row  vector  for  use  in  reduction. 

(5)  Start  loop  to  modify  the  row  vector  by  division  by  the 
pivot  element. 

(6)  Use  inner  loop  to  subtract  the  divided  row  from  all 
other  rows. 
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(7)  When  aLl  row  elements  defined  in  (4)  have  been  modified 
as  in  (5)  go  to  (9),  otherwise  back  to  (5). 

(8)  Increment  counter  of  the  coordinates  to  be  retained. 

(9)  If  aLl  coordinates  have  been  analyzed,  go  to  (10), 
otherwise  back  to  (1)  for  new  coordinate. 

(10)  Use  loop  to  do  reduction  of  triangular  180th  order 
elements  to  be  kept,  to  a  triangular  74th  order  form. 

(11)  Transfer  74th  order  symmetric  form  to  a  74th  order 
square  matrix.  Reason  for  the  double  step  (  (10)  -  (11)  ) 
is  because  of  space  limitations. 

(12)  Return  to  subroutine  STGN. 

4 .  Variable  Arrays 

(*  Signifies  Input  Parameter) 

*  VELCTY(15)  Values  of  reduced  velocities 

A(74,74)  Flexibility  matrix 

0(2,2),  B(2,2)  Transformation  matrices 


Mass  matrix  parameters  for  building 
distributed  mass  matrix  and  also  for 
forming  dumbbell  transformation  matrix 

Mass  matrix  element  values.  Only 
nonzero  values  are  saved. 

Mass  matrix  subscript  value,  respective 
storage  position  stored  columnwise. 

Dynamics  matrix  work  area 

Storage  for  up  to  three  eigenvectors 
for  printout  purposes 

Eigenvalues  computed 

Work  area  for  eigenvector  computation 

Numbers  of  eigenvectors  computed 

Storage  for  variable  eigenvector 
formats 


*  X( 22 ) 

*  Y( 22 ) 

*  M(22 ) 

*  SX( 22 ) 

*  IX( 22 ) 

AM( 150) 

ICT( 150) 

U( 74,74) 
VS( 74 ,3) 

EIG( 74) 
VECT( 74) 
NCB( 3) 
FMT( 3,3) 
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*  i  juWBIWF 


FMA(3) 

*  C0SLMA(4) 

*  BRA( 4 ) 

*  SA( 4 ) 

*  CAPA(4 ) 

*  DELTAY (7,4) 

*  B(  7 , 4) 

*  D( 7 ,4) 

AM(2,2),  BM(2,4)  Work  areas  for  use  in  aero 

CH1(2,4),  PARTA(2,4)  coefficient  evaluation 

Lengths,  bending  stiffnesses,  torsional 
rigidity  for  all  coordinates,  for 
building  stiffness  matrix 

Transformation  angles  and  various 
other  quantities  used  to  set  up 
original  stiffness  matrix 

Subscript  list  of  variables  to  be 
kept  during  reduction  of  180th  order 
stiffness  matrix  to  74th  order 

Reduced  stiffness  matrix 

Same  as  X,Y,M,SX,IX  above 

Transformation  matrix  for  flexibility 
modification 

Upper  triangular  180th  order  stiffness 
matrix 

KMAT(  12 ,12) ,AMAT( 12 ,12 )  Work  areas  for  building  large 
S(12)  stiffness  matrix 

RED(180)  Column  work  area  for  reduction  of  upper 

triangular  stiffness  matrix 

H92921)  Work  area  to  hold  upper  Hessenberg  form 

of  the  dynamics  matrix  during  eigenvalue 
analysis 

MTJLT(  74)  ,SHIFT(  3)  ,  Work  areas  for  eigenvalue- 

INTH( 74) ,NS( 74) ,ING( 74)  eigenvector  analysis 
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*  XDEIB( 5,2) 

*  XDEICC5 ,2) 

*  XDG  J ( 5 , 2 ) 

ALPHA (  30)  ,  PSI  ( 30  ) 
BETA( 30) ,EIB(30) 
EIC(30) ,GJ( 30) , 
AE(  30) ,  LOO) 

IX( 74) 


STIF( 74,74) 

*  PARM(22 ,5) 
FINAL( 74,74) 


REDK( 16290) 


Work  area  for  eigenvector  formats 


Storage  for  surface  variables  input  to 
aero  coefficient  evaluation 


Storage  for  surface  and  strip  variables 
input  to  aero  coefficient  evaluation 


NSZ(74)  Transformation  matrix  for  180th  order 

triple  product 

XTD (3,3), XTE (3,3)  Work  areas  for  above  triple  product 
C.  DFAL18 

1.  Method  of  Solution 


DFAL18  solves  for  the  internal  forces,  moments,  and  torques 
acting  at  a  specific  node  on  a  FPR  wing,  pylon,  blade  system, 
due  to  an  externally  applied  airload.  The  solution  involves 
determining  the  response  of  selected  coordinates  qn(t)  and 
using  the  responses  to  compute  internal  loads. 

The  response  of  each  coordinate  is  calculated  by  adding  the 
time  dependent  transfer  functions  for  each  flexible  mode. 

Each  mode  has  the  unsteady  aerodynamics  associated  with  its 
frequency.  The  transfer  functions  are  expressed  as  polynomials 
in  terms  of  the  Laplacian  operator  s.  While  the  inverse 
Laplacian  could  be  formed  a  finite  differencing  integration 
technique  is  used  to  calculate  the  coordinate  time  responses. 
Internal  loads  are  then  calculated  by  multiplying  the  coordinate 
displacement  times  the  elemental  stiffness  matrix.  This 
product  is  done  at  each  time  step. 

Coordinate  responses  and  internal  loads  are  shown  as  a  function 
of  terms  on  CALCOMP  plots  and,  if  desired,  printed  output. 

The  flow  of  the  computer  program  is  shown  by  a  macroscopic 
flow  diagram  (Figure  15) .  A  brief  description  of  each  sub¬ 
routine  with  a  written  flow  explanation  is  given  on  the 
following  pages.  The  main  variable  arrays  are  defined  at  the 
end  of  this  chapter.  The  reader  is  referred  to  the  FORTRAN 
listing  (Volume  III)  for  program  details. 

2.  Subroutine  Descriptions 

SUBROUTINE  TMF 

This  subroutine  sets  the  plotting  axes,  scales  the  variables, 
and  plots  the  time  gust  function. 

ENTRY  HIST 

This  subroutine  entry  point  prints  and  plots  each  force, 
moment,  or  coordinate  displacement. 
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Figure  15.  DFAL18  Flow  Diagram. 
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SUBROUTINE  RESP 


This  routine  evaluates  the  time  response  of  all  second  order 
transfer  functions  and  accumulates  the  real  part  of  the 
response  for  each  transfer  function  to  form  a  total  forced 
response  for  each  coordinate. 

SUBROUTINE  MASS 


Forms  a  mass  matrix  corresponding  to  the  mass  matrix 
distribution  of  the  generalized  coordinates. 

SUBROUTINE  FN 

This  subroutine  computes  the  time  history  of  the  gust  input 
function.  The  routine  may  be  simply  modified  or  extended  to 
provide  any  gust  input  desired,  according  to  a  user's  needs. 

SUBROUTINE  INVDET 

This  subroutine  computes  the  complex  determinant  value  of  a 
complex  matrix  by  Gaussian  pivotal  reduction. 

SUBROUTINE  MULT 

This  subroutine  forms  a  matrix  product  UA  where  A  is  a  real 
matrix,  in  this  program  a  flex  matrix,  and  U  is  a  complex 
aerodynamic-inertial  matrix. 

SUBROUTINE  DPROD 

This  subroutine  forms  a  matrix  triple  product  ATBA  where  B 
is  a  symmetric  matrix,  that  is,  B  *  BT. 

SUBROUTINE  DFAE ,  BESJ,  BESY 

This  subroutine  and  its  associated  routines  BESJ,  BESY,  form 
the  aerodynamic  influence  coefficient  matrix  corresponding  to 
a  given  reduced  velocity  or  flutter  frequency.  BESJ  computes 
J  Bessel  functions  of  order  0  and  1,  BESY  does  likewise  for 
Y  Bessel  functions.  Methods  used  are  those  given  in  IBM 
SYSTEM/360  SCIENTIFIC  SUBROUTINE  PACKAGE,  VERSION  III 
PROGRAMMER'S  MANUAL ,  H20-205-3,  page  363-364. 
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SUBROUTINE  RESCOM 


This  subroutine  sets  up,  evaluates,  and  reduces  each  dynamics 
matrix  formed.  It  then  evaluates  transfer  function  coefficients 
and  time  responses  for  all  coordinates. 

SUBROUTINE  REDUCE 

This  subroutine  eliminates  undesired  coordinates  by  Gaussian 
reduction,  and  retains  the  modified  matrix  and  forcing  vector 
for  determinant  evaluation. 

SUBROUTINE  ERRSET 

An  IBM  360  subroutine  which  suppressed  underflow  error  print¬ 
outs  when  the  J  Bessel  functions  are  evaluated  very  near  a 
zero  value. 

3.  Flow  Outline 


(1)  Initialize  CALCOMP  plot  buffer  storage  for  plot  output. 

(2)  Start  reading  input  for  a  new  case.  The  termination 
point  for  a  run  when  all  cases  have  been  completed  is 
(23).  Read  all  general  input  and  control  parameters 
and  print. 

(3)  Compute  constants  used  in  iteration  loops. 

(4)  Read  wing  parameters  and  compute  gust  force  magnitude 
for  each  coordinate  station. 

(5)  Read  damping  and  frequency  values  obtained  from  DFAL17 
for  each  reduced  velocity  and  print. 

(6)  Call  subroutine  DFAE  to  initialize,  read  and  print 
aerodynamic  input  parameters. 

(7)  Call  subroutine  MASS  to  read  mass  matrix  inputs  and 
compute  mass  matrix. 

(8)  Read  flexibility  matrix  and  flutter  frequencies  output 
from  DFAL17. 

(9)  Compute  flutter  frequencies  in  X  form  from  input 
frequency  and  damping. 

(10)  Print  flutter  system  roots  and  flexibility  matrix. 

(11)  Convert  flex  matrix  to  foot  units  and  form  modified 
forcing  vector  F* . 
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(12)  Start  iteration  loop  on  all  coordinates  at  which  a  force 
and  bending  moment  or  a  coordinate  response  is  required. 

(13)  Initialize  counters  to  indicate  flutter  frequencies  and 
coordinate  stations  to  be  used  and  call  subroutine  RESGOM. 

(14)  Begin  loop  over  all  time  points.  This  loop  transforms 
all  coordinate  responses  from  the  generalized  form  of  two 
vertical  displacements  to  the  coordinates  representing 
the  dumbbell  model,  that  is,  one  rotational  and  one  ver¬ 
tical  translational  mode  at  each  station. 

(15)  If  only  the  translational-rotational  coordinates  are  to 
be  printed  and  plotted  (input  option  NT),  go  to  (21). 

(16)  Compute  parameters  required  to  compute  forces,  torques, 
moments . 

(17)  Start  iteration  loop  to  evaluate  all  quantities  which 
are  functions  of  the  translational-rotational  coordinate 
responses  computed  in  (14). 

(18)  Compute  bending  slopes  as  functions  of  translations  and 
torsional  angles  as  functions  of  rotations  and  transla¬ 
tions.  If  response  at  the  first  station  is  to  be  evalu¬ 
ated,  set  cantilever  boundary  condition. 

(19)  Compute  forces,  torques,  moments. 

(20)  Call  subroutine  HIST  to  print  and  plot  time  histories  of 
forces  and  moments  computed  in  (19). 

(21)  Call  subroutine  HIST  to  print  and  plot  time  histories  of 
rotational  and  translational  coordinates  computed  in  (14) 
as  well  as  slopes,  etc.,  computed  in  (18),  if  any. 

(22)  End  of  loop  (12).  If  all  stations  of  force  and  displace¬ 
ment  analysis  have  been  completed,  go  to  (2)  for  a  new 
case,  otherwise  to  (13)  for  the  next  station. 

(23)  Close  the  plot  tape  and  stop. 

SUBROUTINE  RESCOM 

(1)  Start  loop  over  all  flutter  frequencies  to  form  reduced 
matrices  and  forcing  functions. 

(2)  Zero  the  dynamics  matrix  and  call  subroutine  DF67  with  a 
reduced  velocity  given  by  the  forward  speed  and  the 
flutter  frequency.  Upon  return,  the  complex  aerodynamics 
matrix  transpose  is  provided. 
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(3)  Add  the  mass  matrix  transpose  to  the  aerodynamics  matrix 
transpose  to  form  an  effective  mass  matrix. 

(4)  Gall  subroutine  MULT  to  form  the  product  of  the  flexibil¬ 
ity  transpose  and  effective  mass  matrix.  The  result  is 
the  transpose  of  the  system  dynamics  matrix. 

(5)  Transpose  the  dynamics  matrix  to  get  in  normal  column 
form. 

(6)  Choose  the  particular  flutter  frequency  defined  by  loop 
Cl)  and  substitute  the  value  for  X  in  the  matrix  A  -  XI. 

In  addition,  form  the  complex  valued  forcing  function. 

(7)  Gall  subroutine  REDUCE  to  eliminate  all  unwanted  coordi¬ 
nates  and  retain  the  6x6  matrix  and  6x1  force  vector 
necessary  to  compute  determinant  values.  If  all  coordi¬ 
nates  have  not  been  done,  go  back  to  (1). 

(8)  Call  subroutine  FN  to  build  the  form  of  the  time  gust 
function. 

(9)  Start  loop  over  all  coordinates  to  compute  transfer 
function  coefficients  for  each  flutter  frequency  and 
coordinate. 

(10)  Pick  up  the  reduced  matrix  corresponding  to  a  particular 
flutter  frequency  and  substitute  the  forcing  column  for 
the  coordinate  whose  response  is  desired. 

(11)  Call  subroutine  INVDET  to  evaluate  determinant  of 
substituted  matrix. 

(12)  Store  coefficients  corresponding  to  a  particular  coordinate 
and  frequency.  If  more  frequencies  remain  to  be  evaluated, 
return  to  (10)  for  another  flutter  frequency. 

(13)  Call  subroutine  RESP  to  compute  the  response  of  a  partic¬ 
ular  coordinate  as  a  linear  combination  of  the  responses 
of  each  flutter  frequency  transfer  function  to  the  gust 
input.  If  all  coordinate  responses  have  not  been 
computed,  return  to  (9)  for  another  coordinate. 

REDUCE 

(1)  Set  up  column  subscript  vector. 

(2)  Use  loop  to  divide  determinant  by  denominator  factors 
formed  when  the  first  six  flutter  frequency  terms  are 
evaluated  at  a  given  flutter  frequency. 
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(3)  Start  loop  over  all  coordinates  to  be  reduced  or 
eliminated . 

(4)  Use  loop  to  find  the  largest  diagonal  or  pivot  element 
remaining  of  those  coordinates  which  are  to  be  reduced. 

(5)  If  all  coordinates  to  be  reduced  have  been  reduced,  go 
to  (14). 

(6)  -Set  indicator  for  the  coordinate  to  be  reduced  to  zero, 

multiply  determinant  by  the  pivot  element  chosen  and 
divide  it  by  the  corresponding  denominator  element 
evaluated  at  the  given  flutter  frequency. 

(7)  Use  loop  to  divide  forcing  vector  and  matrix  row  by  the 
pivot  element. 

(8)  Set  the  reduced  force  element  to  zero  to  minimize  compu¬ 
tational  error. 

(9)  Start  loop  to  reduce  matrix  coordinate. 

(10)  Skip  to  (12)  if  product  element  is  zero  or  row  is  the  row 
to  be  reduced. 

(11)  Use  inner  loop  to  subtract  the  modified  divided  row  from 
each  row  defined  by  (9). 

(12)  If  all  rows  have  been  modified  go  to  (13),  otherwise 
back  to  (9)  for  new  row. 

(13)  Use  loop  to  zero  the  reduced  row  and  column  elements. 

Go  to  (5)  for  new  coordinate. 

(14)  Start  loop  to  pick  out  those  elements  to  be  retained 
for  the  6x6  matrix. 

(15)  Pick  out  all  elements  of  force  vector  and  matrix  to  be 
saved  and  save  for  column  substitution. 

(16)  Return  to  subroutine  RESCOM. 

SUBROUTINE  FN 

(1)  Set  the  time  zero  forcing  value  to  zero. 

(2)  If  the  gust  number  is  1  go  to  (4),  if  greater  than  2  go 
to  (7). 

(3)  Compute  number  of  points  for  half  cycle  sine  wave.  Gust 
form  is  sine  wave  with  input  frequency. 
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(4/  Compute  sine  wave  using  recursive  re lations hip  for  sine: 
sin('HAiJf)  =  siny  cosA ip  +  cos^  sinAy 

l5)  If  the  gust  number  is  1,  go  to  (II). 

(b)  Use  a  loop  to  square  the  sine  wave  for  gust  number  2. 

Then  go  to  (.10) . 

(7)  If  the.  gust  number  is  not  equal  to  5,  go  to  (9). 

(8)  Compare  Kussner  gust  function  as  function  of  time,  forward 
velocity,  and  semichord.  Then  go  to  (11). 

(9)  For  gust  numbers  3  and  4  a  step  input  must  be  computed. 

Set  the  time  function  to  1.  If  gust  number  is  3,  go  to 

(11). 


(10)  For  pulse  gust  and  half  cycle  sine2  wave  the  remainder 
of  the  response  must  be  set  to  zero. 

(11)  Call  subroutine  TMF  to  initialize  and  plot  forcing 
function  time  response. 

(12)  Return  to  subroutine  RESCOM. 

4 .  Variable  Arrays 

(*  signifies  input  parameter) 

Complex  natural  frequencies  in  S  domain 
Wind  force  vector  amplitudes 
Coordinate  response  as  f(time) 

Flexibility  matrix 
Flutter  damping  coefficients 
Modified  gust  forcing  vector  amplitudes 
Flutter  frequency  coefficients 
Nonzero  mass  matrix  values 
Storage  for  transfer  function  coefficients 
Storage  for  coordinate  time  responses 


PC( 10 , 2 ) 

FV ( 7  4  ) 

Y  (1003) 

*  AWK(74 , 74  ) 

*  G( 10) 

FY ( 7  4  ) 

*  0(10) 

AM( 150) 

XM( 10 , 6 , 4 ) 
YC(1003,6) 
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G  (2,2,22) 

Trans format ion  matrices  for  each  linear  to 
displacement  rotational  coordinate  transform 

MSB(150) 

Nonzero  mass  matrix  indices 

NSG(22) 

Starting  coordinate  numbers 

NSB(4 ) 

Numbers  of  coordinate  response  to  be  run 

NSGB(10) 

Coordinate  numbers  for  which  response  is 
desired 

BUFC512) 

Buffer  storage  for  plot  output,  multiple 
of  256 

TTLC20) 

Run  title 

DY( 22 ) 

Section  width 

BC  22  ) 

Section  length 

XDC1003, 3) 

Storage  for  responses  to  be  plotted 

El (10) 

EIB  section  values 

GK(IO) 

GJ  section  values 

SLG(IO) 

LI  section  length  values 

SLD(IO) 

L2  section  length  values 

A(74 , 74  ) 

Complex  flutter  equations  matrix 

Z(74 ) 

Work  area 

*  Q(74) 

Zero  velocity  poles  or  flexible  frequencies 
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APPENDIX  I 


SAMPLE  PROBLEMS 


Sample  problems  are  provided  to  assist  the  user  in  checking 
out  and  using  the  computer  programs.  Each  sample  problem 
consists  of  the  data  sheets,  the  input  display  included  in 
the  output,  and  condensed  output.  Examples  of  the  CALCOMP 
machine  plots  for  programs  ARAP06  and  DFAL18  are  also 
included.  The  sample  problems  are  taken  from  the  D270A  FPR 
dynamic  analysis,  discussed  in  Section  IV  of  Volume  I. 


FORTRAN  Cm*  i*g  Farm 
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ARAP06  Sample  Problem,  Input  Data 


A RAP 06  Sample  Problem,  Input  Data  (Continued) 
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ARAP06  Sample  Problenv  Input  Data  (Continued) 


■**  O  it  —  vr»  ro  ^  o-',l 

rsj  •£>  OQ  v£  o  q«*»o- 


—  w 

.tf  v 

i  -  j>~ 

->  W-  V/  <*i  <5- 

^  M 

o  o 

-  4  r< 

—  —  o-  *> 1  •*' 

r*  r>4 

• 

*•  -  •  . 

—  »  )  - O 

1 

r 

i  i 

i  i  > 

« 

r—  ^ 

**  "4- 

lOcrTVl  .r> 

0“  ^ 

-*  M 

,~-ct 

o  <r  g-  o-  ^ 

o  o 

M  o  Vrt 

—  o  r»  O  O 

j  »'  ' 

1 

i  i  i 

ii  i 

>d  Vj 

<0  V»  <» 

M  v»  r-(r  ••  x 

95 


ARAP06  Sample  Problem,  Input  Dat  (Continued) 


FORTRAN  Coding  Form 
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ARAP06  Sample  Problem,  Input  Data  (Continued) 


FORTRAN  CWmg  Farm 


98 


vs  vH-ol  ' 

*  "T-+ ;£<■-  o 

'  €*>  Uj  *1  A,  v»  Q 


<n  -  ^a(  " 

IN,  JO->1 


-*-  V)I«M^ 

;S5-^jK« 

^tWs'T 

°  »"  °®  o  *'*!'» 

;_i_ulil._ 

•  bJo-1—  c\A  M 

^  '**.  vi  O 

I  . 


i«jrl-  «~!  |V 

-I  A  H  f-~  r"j  r< 

*!_§'  tL^I—  .  '“*>C 

5  »!  '  ’  l~  1  '  'I 

I  '  -J- o  vS  *o  r'  W 

3 

I  ;  ,  o,  v-  ^  »  n  ~a 

1  uiLJii 

r»  J  —-*- ^ 

'OqOm^;  o  -,-~h 

:  .  1 11 1 1 1 

^i5o  r  n  w 

>i0D  i)  o  m 

■  ]•■•!•-;  •  •' 

Q  «  Q  i*»1fll'’ 

V> 

1  1  i  j 

"i,  6q  Q  *S  -^>,  O  ~~* 


snsli  s;^ 

,°°  JO,  <B  <?2  O 

.sxm 


X  Vit  V,  £  vJlJVo 

/O  ift  a  .  .  v\  ■  <r  • 


b#  r-  sS  N  i"  V)  i*1  bo  **)  ^-  (-~1 0® 

>0  ^  NJ(  K  >^nS 

>-0o'o'ONi  <3  V'JsV)  n  n3>  r- 
<A  N  !h<>  n  CltoO-miNJ-^. 

fnN  ^  |  ft)  | 


TiPTmlTR- 

io  S  ■*■  rtV  fco  J)  If)  N ,  i*1 
■'•'vO  >\s'  Oa  i£  *o  fv'^ 

f'sS'i-  O  m  N  II)  r(  HI  .  If)j 


>-' V^-iNi*  O*.*  --■  wTO-'K-’V) 

K  I*}  1*1  ~  ‘•o  '*  r*:> 

bo  <r>  in.  ^  >  H  r, 

js  N  «n  fO  ^  q  V)  o  N  m  n- 
Jv^tvO^ , in  I*), lA.Vl^ >|h.« 


-:'t  M  I! 

:  j  :l  a,:  j  . 


>  h- V)  r-  O  Q  V)  if)  m  6c  rg\j  : 
^  o  d  \S  >n  ">  N»  O'  Vl '  n  I  N  m 
(0*«O  HV  'Co"  jn'ug'v9‘«-  (*> 

<j~  <0  i*i  n  i^i  in'  r-'  ni  m;  bj 

>S  bo  m  i*ir)N  m  |  N  i  if) 

I?!  1  I  I  I  I  I!  I  ■ 


-»  Oo  bo  o'  1-  On' m1  (nJ  o  -j,1  *11 

tv  nrj  V)  M.tfcNs 

rfl  o  O^oo  N  r'i’^lfr'vS 

.  isLst.^N^ ;  J  JJSJW  Jras  <ai  f-i 

qt-ns  m  m  r(|in>  m,  v> 

>'  ,  -  I  t  III  ! 

•ti  ,1  J  J  .:  ,.  J  ,1  j  j  J  , 

bo  |Q  ^  bo 

<J >>  N  r,  r,  r'j  ^'lo  N  q! 
-3-  6o  bo  —  o  a  «i  n  ^'4  n  1 
■*->  9so  in"^>  m  <fv  o  Ci  -  s° 

.  1*1.  1  .  i  1 V  i  _i _ 


V)rO'obOvaV)*Oir)(f)NboO 
V  vS  bo  !><■  'S  <N  ‘*1  f'  f*  f>0  r~  1*1 
bo^<*lQvS's±.UcGON&]*air) 
.S  rj’ao^*  rvO*> 

— ■  ^  •  .V_  K.'  lA  M  ^  I  .  ...  I 


<r^<o  o|^> 

:r*  K  , 


U  nl  bo!<V'^, 

_ t-_ 

-|IN  vL:  Q  — j  I/)  a 


m  m,  (\ 


l  •!  .'  .•  .•  .'  .'  .  ,  .  .1  .  . 

to  d  -V^  m’Vj ; p~  n 

n  rs  r^m  «s  >*)  "^ 

r-  IN'}-  *M  IN  '->->-  Q  »0  >  o 

oVi!~->'mjOi>vjt)'tr-foV)ir' 

>a  v>;fb  O',  i  ■  cm,  m,  ri  -i  hti 


t  at 

T  a*  £ 

^  h 

^  Ul  Vb 

sit 


151 


-*  n  (O^  — ,  bo1  -3-i  r~l  Oo  _  »• 

VTOsjsaas  - 

n  h»  h  (y^d>c  %S  bo  b*  in  *  i 
Sv  -» 1*4-  —  CQ^  p»  ”  boNJ  o  ^  • ! 

vs:  H  o,  i  mi  if\  -H  i-  :  i 


99 


ARAP06  Sample  Problem,  Input  Data  (Continued) 
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ARAP06  Sample  Problem,  Input  Data  (Continued) 


JUNE  II,  1971  SAMPLE  CASE  M*  PROGRAM  ARAP06  USING  MODEL  0270  PARAMETERS 
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ARAP06  Sample  Problem,  Input  Display 
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ARAP06  Sample  Problem,  Input  Display  (Continued) 
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ARAP06  Sample  Problem,  Input  Display  (Continued) 
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ARAP06  Sample  Problem,  Input  Display  (Continued) 
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,  Input  Display  (Continued) 
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ARAP06  Sample  Problem,  Input  Display  (Continued) 
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ARAP06  Sample  Problem,  Input  Display  (Continued) 
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ARAP06  Sample  Problem,  Input  Display  (Continued) 
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ARAP06  Sample  Problem,  Input  Display  (Continued) 
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ARAP06  Sample  Problem,  Condensed  Output 
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ARAP06  Sample  Problem,  Condensed  Output  (Continued) 
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ARAP06  Sample  Problem,  Condensed  Output  (Continued) 
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ARAP06  Sample  Problem,  Condensed  Output  (Concluded) 
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ARAP06  Sample  Problem,  CALCOMP  Plots  (Continued) 
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DFAL17  Sample  Problem,  Input  Data 
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DFAL17  Sample  Problem,  Input  Data  (Continued) 


DFAL17  Sample  Problem,  Input  Data  (Concluded) 
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DFAL17  Sample  Problem,  Input  Display 
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DFAL17  Sample  Problem,  Input  Display  (Continued) 
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DFA1.17  Sample  Problem,  Inprt  Display  (Concluded;) 
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DFAI.17  Sample  Problem,  Condensed  Outpu 
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DFAL17  Sample  Problem,  Condensed  Output  (Continued) 
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DFAL17  Sample  Problem,  Condensed  Output  (Concluded) 
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DFALI h  Sample  Problem,  Input  Data 
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DFAL18  Sample  Problem,  Input  Data  (Continued) 
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DFAL18  Sample  Problem,  Input  Data  (Continued) 
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DFAL18  Sample  Problem,  Input  Data  (Concluded) 


Oils)  I  CAife  ilNc-Syu*Kt^  ic  M  *0  futo  17*  **0?S 


3 


o  onmonoQ^onrioooiior  '*»<■« 


✓>  KJttMUKuiwaiwaiwM.thtntiXJwwMMwit.iu^ 

o  ooaoooooocioooo«>ojoooov 
a  o^juoooouoooooooc»oor*f>i>o 

2«r*  fv  a  ®  (*  ^  o  *m  o  o  o  «*  o  d  o  *r<  ^  o  o  >j 
r~  -t  xinrouiri^uoHioy  t>  o  ^  v  u  O 

O  -  —  —  *»  "M  '« 

*  o  o  o  o  o  o  r*  Li  o  i j  o  o  o  •->  o  n  '  >  n  o  r»  o 


DFAL18  Sample  Problem,  Input  Display 
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DFAL18  Sample  Problem,  Input  Display  (Continued) 
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7  ©  f«>  -,  •»  7  N  «  u  It  it  r-  s  *"  ^  S  O  —  p  «s  *•  lit  N  -  POP  «/  O  ^  O  N  N  /'  #  O  sP<3  — 
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DFAL18  Sample  Problem,  Input  Display  (Continued) 
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DFAL18  Sample  Problem,  Input  Display  (Concluded) 
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DFAL18  Sample  Problem,  Condensed  Output 
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DFAL18  Sample  Probler  CALCOMP  Plots 
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DFAL18  Sample  Problem,  CALCOMP  Plots  (Concluded) 
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ARAP06  FORTRAN  Listings 
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ARAP06  FORTRAN  Listings  (Continued 
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ARAP06  FORTRAN  Listings 
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ABAP06  FORTRAN  Listings  (Continued) 
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—  —  O  uftu  —  —  —  X  **tv  —  o  v  •  ►  hi  o  —  —  •  -  m 

Xft*  44—  h  x  X  —  ft  —  •  •  ft  X  —  Oh  V  m  • r  o  *h* 

4  4—  O  ft  in  x  —  —  v  X  4  V  h  -  ft  a  •  1  Vh*h  4  *  ft  * 

•  ft  •  •  4  INSOlTiaa  ft*—  ft  C  ft  ft  •  ft  •  «C  O  om 

—  —  •  h  —  *  x  X  X  ft  »  — —  tuC4*w*rha*g»  •  h  4 

ft «.  o  saw  •*•  —  4  >«  ««  on  *•  •—  —  him  •  *  —  h  u 

•  s—  c  ft  o  —  —  —  —  ,  ft  ft  04  **•«••  —  — «  •—as 

C<;  V  •  ft  U  ft  ■#  •  ft  •  VI  M  I  UN'  •  a*N  N  •  —  ft  a 

—  C  ft  —  —  —  —  4  C  4  h  ■  •  *  —  h  »  U  a*  X  • 

—  ft  4  -  ui  —  n  n  n  \Mh  -  x  «  —  a  hA  he  —  h< 

•M  —  ft  ft  X  *4**»  *  —  h  40  *  4  V  •— O  •»  4  VI  4  • 

XX  4  X— --mv*XC4  4  4  •  ft  r  •  h  h  ft  ft  IS*  •  ft4ft 

£  ft  ft  —  —  —  X  —  m  *  — «  4  *  —  i*  —  —  h  ft  —  04 

•  X*  Xw»4h  e  •  hft-i*  4  •**h*  •*►_.• 

—  Oh  •  X  X  X  —  —  tnfthftvO***  •  4  o  a 

ft  •  O  —  44  4  Mh  ft  —  —  -—ha*  ft  •  »  •  ft  a  -a  X  ■ 

—  —  •  ft  ft—  X  M  0*hh  —  a  sCCN  ft  ft  «  4 

V  04  —XXX  U ►  ft  ft  X ®  •  •  V  —  4  —  4  X  4  4  O 

4  V  —  ►  whhft)  ft  4  »  •  —  —  ft  •  •  *  fw  ►  *-  •  ft  U  •  ft  •  ft 

•  —  •  »n  COlh  •  4  •  •  •  «  *»  X  •  » 

—  hv—  Xvw  a  Maa  —  —  •  •  Wl*ftft**ac*h«*^e 

ft  _ «  5  o  4***  m  •  c  •  ftftfto**vh*e«<N(v  x  O  v 

—  e  i  a  —  xvw  4  w  —  —  —  ft  —  —  —  v  v  k  •  c  —  —  *«*  4  •  e  v  v  • 

►  m  C  •  —  X  •  h  ft  —  4  ftp04vrv  •  a  »  »  -  •  O  4  v  v 

—  4  w  m  —  ft  O  HV—V4V*4  —  ¥■•  «*  V  —  •  >  a  ►  vt  V 

ft  aUSr  Oaa  XXXhrVVX  ---aa--  a  al)  X  •  •  —  V  *  -  h  -  - 

—  —  •—  *ftl444Cftft  fth**uO  •  V  •  • 

ft  ft  ft  X  —  XX  XXXXft  •  ft  S  ►►►►►*  -  •  *4  u4f  *  —  hah  ► 

4  —  4  C  4  ft  Ui  •  •  »  4  h.  h  ft  ftftftftftftfthOh  hh®4*ftftft  ft 

X  OhlNJl  ft  S  N  S-  -I  -  X  t  *  X  I  X  ►  4  »  O-— IM-ftfftt  X 

4  4  •  &  a  ft  —  h***44*ft  44ft4449*h*X**h*44hi  « 

C  O  ft  —  —  —  4  ft  V  UUOOOOft  ft  COO 

U  O  U*  4  4  4  4  4  4  X  444 

—  hahft  —  fM  —  -Nm-NN*  —  ha  ft  ft  ft  «  h  •  ft  H  —  h* 

OM  ft  «n  ftp  N  f  N 

u  —  —  —  —  —  ft  ft  ft  ft 
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FtU 


C  V. 


4  ♦ 


•  « 

iioS.* 

4  w  b-  u 

•  Q  ui  * 

•  * 

•  •  •  V 

C»  •  •  * 
«Nll» 

►*  IUIU 

4  a 

•  u.  V  • 

•  •  1/1  • 
*r>  -  a.  rv 
ft  •  w  • 
u  4 

U  »  ft 

•  -#  •  • 

•  «  Hi  m 
N*  •  I 

•  *< 

«  «  U.  </> 
IbN  *U 

•  *"  •  a 

;  .'g 

■  * 

**  o  o  • 


o  • 

u.  • 

c  —  - 


*  <r  •  • 

•  c  ■#  • 

>  —  a  *  i 


•  •  *  •  M 

al  •  N  • 
e-*  • 
n  ui  •  a 
•  0*0 


s .  : 


£5 

V*  4  - 


tr  b»  #  X  • 
ft  4  •  U 
•  »  U.  •  Uj 
•  m  •  at/’ 


•n  ft  4 
*-  o  ••* 


2*y.** 

►»  7  •  Ui  • 

•  >QN 

rs.^ 

*»  h,  Hi  • 

•  •  • 


•  ••••40u. 
ONOVUblk 

a  *  m.  ui  uj  «  o  •  ft 
►  ■  O 

•  U.  ft  N  •  • 

•  •  *  >/)  *  m  •  x<  • 

*  •  •  ft  ft  o  m  a  c 

a  «  o  o  O  *-  •  a 

t_>  u>  O  H-  •  ft 

i/i  C  •  •  «  lb 


•  Fli 

•  Ui 

«  r 


•  K  «  •  Ui  • 

a£2o?5 

-  !-i.  :->i 


©  — ft  ft  * 


4  a-  « a  • 
• «  HI 
4  •  4  ft  •  • 


■  ••ki*  UH  atfiyui  •  • 

• 

MO«D  • a  • 4  * %  • 4 
•  bxa  •B?bix\bft 

Z  \  • 

XV  a 

UiaM^Vf  a.tfCVb  •  Ui  — 

X  ft  Ui  <  ft  x*a  • 

ft  Ub*  a  ft  •  «i  •  H  4 

a  •  ft  •  b  .  «  •  xa  C  •  ft 

y 

• 

4  |  Ui  *x  Hi  b  *arv 
ft  e Mft  -  ••  ►  ft  m 
»li*  *  “  •  — 

•  >  •  •  a  X  *  aa 

Pi 

• 

•  • 

hi 

•  •  *u  *  *4  ®  *  b>  •  •  ft  •  • 

*  *C  4  K  Ht  •  ft  *•  4  •  Hi  •  x 

Ift  •  «  • 

V4 •  xbnex 

■^4 

•XU 

"blx 

•  ONC  •, 


•  •  m  •  4  j 

4  •  it  u.  pOil 

u.  ft  4  •  •  •  JUJ  •  < 
•NIL  •  •  W  O  • 

•  •  i  i  v  r  ■  ' 

■  •  ■  s«  ■  » 

d i •  •  « 
uw  a  a  4  Hi  ft  « 

C  Cl  i/>  <  Cl  ►  r  •  Ui  a 

t  -t  r 

O  •  a  O  a  >  ^  u.  a  ^ 


Hi  4  tn 


v  •  •  •  a  •  u.  • 

WMI*  m  •  •  4 

•  «o  i  a  ••  •  it 


*  k  o 
.  **.; 


•  •  o  • 
Z  ui  • 

HlbO  I 


•  ft  •  •  i  n  •**  jUirc  r  •  *- «  *ft 
9*«a  «•  ►  «ZK"fi*  u  •  ft  • 

*  —  a  •  *  *  •  — 

.  _  _ _  »if>»  rtiOMi 

•  Ct.  C  M  1*  I  *1  •  4  a  Ui  ft  IX  •  ft  —  U  •  •»  «- 

•  •iUjMXEMt  a  •  •  Hi  •  a  a  ^  •  -»  a  •  ft  * 

I  h  Vi  I  lb  ft  4  litl  ;ft  V  •  •  Hi  ft  —  Z  -t*r  *1 

ft  Ob  i«l  —  -  - -  --- 


_  M  «N  ••  S4  »  ><\- 

w*  *o»  *■  a  •  o  s.  v  ftft»ftftHvft«ftfteft*» 

“  I  KHSMM.il  ^  »\|  •  ■  v  XT  T  X<II\ 

r  a  w  a *•  U  «* •  •*  •  •*  x  •  •  -  x *--• 

•  m  •  m  m  IS  t/»  |  •  UI  •  •  ft  a  •  •  S 

-  -  *•  C  ft 


>  a  v  s  •  4 

•  *  •  «•  o  «»  *  O  •  a  w  •  *  y  » 

Ux^  *  •  *W  *(T  af. 

Ui*  aH  ab-Na  NtbxibXM  »  •  t-  »  b-  •  b  N  ►  UI  b  Z  X  ft  ►  *-  *  aC  b  •  *  •  C  b- 
O  aNii<  3  *04  •  V  Ui  i4aC>44444  i«0«»  »  a.lb  «  «  NF  «  «  N  <  f  4  « 

",Icscl2'J* I5,  '“"ifj!8ts,8s:‘,,,!r“gt'J  e 

r  »  ^  ^  xv  fu  a  a  ^  4  r»  a  r  k*«*  iu  ^  y‘au*fc»«,i,Hi,fctu»ua4,fc,fc#»Hiib,k‘-*^H.a,fc 


S  ft 


A 

Q  (v 


5 

I 

O  C 


*  s 


«  M  a* 
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l*»  fWMl  «/•  KIE.III  ICMUC- F|0.«.tMI.«l 
14#  FORMAT  (/•  mill  I  INI'.  II.  10FU.21 


r 


O  u 


0  O  G 


i 


•  (A 

t~  Mi  • 

►  V>  ► 

-  3  S 

«  T  • 

►  •  X 

-  O  4. 

i  7  «. 

:  i  : 

•  •  Mi 

«  U  M 

£  iS 

c  52  *„ 

'  fK  e 


s  I 


=  2“ 


•  ■  PM  P* 


.  sr 

►  x  >  *  • 

1*i  ?  w  *  £ 


5  jj  B  '« 

•  a  •  Mi 

a  3  a  i 

•  •  «  v* 

<  *4  •  W 

«  *  ►■VI  «* 

PM  ►  VI  <  PM  • 

36  i*t  s 

« ►  H  »U  « 

•  VI  O  J  —  » 

r;  am  ir 

4.  4.  •  a.  VI  c 

4  O  p*  «  «?  • 

I^O^Wf  Mi 

in  vi  mi  Vi  «  C 
*  •  >  o  •  • 

4  <f  V 

«s  e  v>  u  •  C 

r  •  *  •  w  • 

4.  V'  ►  X  J  v> 


c.  .  ..  X  • 

r--f  **  o 
u  •  O  • 


jg:sjs 

r  *k  • 

.  *M  O  >  •  U 

-IU4 

s-sis- 


V  -  V 

if*. 

*  ip 

S«S 


Eyj 


©  J  1  • 

S*»?f s»ss * 

£_-5j_-45jI-^e 

m  I  a  I  I  ■  I  I  I  ■  I  f  e 

•■“••••■••■•■••••■••■PR*-* 
!•>•••••••••  •  © 


■  -  «|M  I 

•  ►  NOI 

o  —  r  ►  ► 

,4.  *4.  VI  • 
*-**«  •* 

•  pm  •  M  «  | 

VI  ►  4  X 

•  •>  a  o  o  I 
►  ••.►x 

DNKV I  ► 


l-'^S 
\zZSZX 


7 

* 

>  ? 
v«e 


♦  X  VI 


o* 

K 

3  ! 


4  Ml  •  PM 

*  8  1  = 


«U  •  U 
f  4*  • 

lit  J 


t 

tZu.i  K 

m  mi  «  •  •  —  ► 

■°i;s 


~  I*  U 


I  2 


•  »IUKMIp4« 
Mi  O  ►-  4  ►  V  • 

►►«4»«©^4 


mw«.p«wv.p*at 

>^M«*UWWUIM»00  _ _ _ 

•  I*rr-a»4  «.  «r  •  i>  u-pypf  i)4V 

S  •  H *  *  «•  »  •  u  S 1  •  ►  **  r  ••  X -M  * 

"  O O  »OOM« 


nt-PO  rOa»Hir<\irrBBiMiW'»'inim'»4ir  •  O 


iO*OOf V 


Vi  - 

•  o 


r=-5T-5 

O  ►  O  —  Cl  Ml 
«"« 

Mi  Ml  M  U  Mi 

m  m  —  x  « 


■  mh-pppM-mppIui  m  •  C  Ml  M  Ml  Mi  N  •  •  •  Mi  Ml  lb  C  Mi 

•  •••••  •  «  I  •  »  »•>  I*I««4V  »r  •  •!«««•«  •  •  •  *  • 

nr  aMMnriMrriNtirr  <  im>4ip4»i>«iik«vi'‘ 


JmOOOCI 


D-rWluhMppp-lMCwOVPP-a-Kiiiip 
COOCUOaoCUOMM«MK^^Mp«wMP<P^^^b'b>>PMPM.B»Bv-4 
«4«4««44444  »  —  p*  -■  «  ■■  mm  O  C«I 

*  «  4  4  4  4  4  4  *  4  4  4  5  5  *•  — “  •-  lSlSl  -  —  *  —  5 


IE 


*  • 

• «» 
►  •• 

fw 

•  4 

2- 

•  PM 
►  •• 

l: 


X  ~ 

r. 


*r 

►■  x 


s: 

*  *-  •  D  — 

«  P  Mi  Mi  V  Ml  •  4 

•  UV  C  • 

*  «  —  p.  C  PM  V  — 
4  •  m  ■  s  m  •  » 

iniioi  vr  «p 


<ooor  n  4  ■» 

C  Mi  »w 

•  •  4  » 

vi  •  »  e 

u  ►  v>  p 


T3 

Q) 

3 

C 

•H 

-P 

c 

o 

u 


m 

Cn 

c 

•H 

■p 

to 


H 

« 

g 

fci 


VO  I 

O  i 


5  3  3* 


A 

O  c 


p 


»4  o  1 
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▼ 


t 


o  o  c 


▼ 


e 

JS? 


-  -  a 

l  •  W* 

-  -  * 


IW  X  •  T  •  ft 


.i 

~  wi 

•  X 
■ft  • 

•  mm 


5557 


X  2 

ft.  m 


51 

s  % 

-•5 

•  X 


■XXX*  - < 


a 

t 

* 

? 

o 

G 

-  » 

5  s:‘ 


1  c  • 

_•  It 


52  'I  2 


ft  *> 

5-* 

:-s 

•  X 


—  ft  ~  - 

•  N  t  I 

«  *iz 

•  X  X  «ft 


ft  ft  14  "» 
•  2  • 
—  a  — 
*  •  n 


x  ^ 


ft  ft  -»•  • 

eS;::::-;; 


►  X  -ft 
ft  **  •  ft 

..  •  •  X  -  X 
ft-  —  ~nnft«e-»  •  ft  »/> 
t  -  a  -  (X^Oulvi  •« 


til 

«fl  *  Vt 


?3S 

XXX 


s. 

■?* 


•  *X  ■ 


^ 

•  •  •  •  m  I 


X  *  -« 
•  IT  ft  ft 
••  p*  X  ■ 

iiS. 


f 


«  X  X  x -- - 
M  ft  ft  ft  V  ft  ■» 


sn; 


;ss 


—  n 

tu. 


•  •••  •  «k,  •  •  o  uftftiihftftftzxzoftoir  •  r  •  t  •  ► 

«r  *  *  irif— rf"xn-»c  ftc«4ftroiOftO» 

*w.r~iCrXXXXXX-".-«'X--M~ft  »ft  •  ft  X  1 


or  •*•( 


soSSSS  «o w 

«di tftftKC* 
—  _  X  X  X  X  —  £  ft 

*  M  •  ftftftftft*'* 

BP>OXXrXft  ft*r 
■  »  ft  X  K  r> 

Mir  «•< 


=  ^  = . i-!£- 


_  _  _ _ ftftftft  —  ft«-XHI 

•  +  t  «  —  —  XXXXXX-»»-  —  —  —  —  ••—  ft  ft 


ft.  *  9  <•  W  -  X 


ftt  ft-  -  X  X  X  X  X  X —  ft  ft  »• 

u.ftwc«*.ftftftftftftftftcft«ftft«««ft  ft  ft  ft  I 

•.  x  •  a  -  a  x  2  x  x  x  x  a  *  a  a  a  *  a  a ~  -  »f  ' 


sltut;*:**! 

*  —  *  T.T  Z  X  >•  . 


*  it 


u 
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MCE 
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CT  c  c  c 


O  O  O  t_ 


I  i  ;  ! 

l'|i 


2i  -  « 

32  2 1 

•  •  ik  • 


uc 

*•  —  -►  I 

•  —  4.  #■  3 

Cl  *  -c 

•  ►  w  *  • 

•  Sue 
» •»  • 

2  m  —  -o 
*  *  «  o 
« c  /►/ 
i4m  •>*  tr>  • 
ft -  «•«*  I 

3?  ©-It  *  * 

5v.8.5.i 

• 

•  •  «  v  4  c  r. 

X  »•  •*»  -  «  X 
OMOIOf  • 

U 

«A  •  «.  4  W1  S  X  fc 

ft  c  -  *ftZ  •  5 

•  X  *•  —  ft  ft  •  X 
X  *-  •  «  «  u.  o  • 

►  *  •  -  -  ©  —  %* 
~3CftU+«ft 
O  -  »  I  4  2  * 

fgf.*:-* 

S  *  v*  -  w  « 5 

►  "7^Z  •  *  • 

M  MN»*<  M 

£  •  e  c  <r  f  *  o 


s:U£:?| 

2 

5wc  c  u 


=»  5 

5825 

2522 


"irrssw*  £! 

14-  « 


M  k  *4  m  m  m 

iff?  tt t 

— 25 

!  <*  «  ©  14  Ik  (4 

,2SS  ,  .  , 
SSS  -2- 

1  14  14  14,  OOP 


I-;;  kvk 
eeo  •*  —  » 

ft  —  X 
♦  ♦  ♦  •  w  ft 

«  «  «  14  14  14 

W  *  *  -  »  - 

ft— X  ♦  ♦  ♦ 

OCH4 

©  <4  14  «4  —  •  “ 


8  —  ft  ft.  ft  •  14 

•  I  I  *  ft 

*4  •  u  ...  C 

*y*X©—ftU  u  ft  I  O-C 
ft  Z4«It  I  •yuZlOW«lA> 
•-  <<*«5l4t-3i«ftC 

1  I  ICt-»  ©  ft  -  *  1 

e  o  c  —  —  t  •" 

*i¥**¥¥eu&>a  *  Cft  o  *  •»  u 
ft4ftft44C«vlllOfi>*UXO« 


»  1  ♦  * 
*r  S-  3 
£* 

CMItfft  - 

9ef 1 aJ e 
1  ©  •  ■  —  © 

)0«J4CC 

IOV4  xr o 


*1  ft-  -  x  I  • 

ft  —  —  ©  *  ft  ! 

♦  ««<Oft-ft.ft  K 
KKtr«»>» 

rs=x,...  ♦ 

©  ©  ft  ►  -ft  « 

«44ik«*>r  3  • 

m  u  r  z  •  r  *  m 

>  1  •  • 

■  1 3 • 

1*4  -  -  I  Z  4 

«ft-iftftft-XOftU 
^►V4ft*ii4Uftft 


se  s; 

—  —i  14  V 


ii;s3  .. 

TT5i i?I *4 
*  «j*  *  u  ©  •  »5 
—  «  w  ft  t  ar 

eo  ••  ••-ft 

'••.w’tf S2 

?8:2?:tsbr 


a:  i 
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CARO  COUNT 


ARAPO 6_ FORTRAN  Listings  (Coirtinued) 


▼ 


w,  y  U  Q  ^ 


<r  «  inr  ► 

IN  M  m  m  *■  •  • 

•  •  •  lO  *\  • 

m  FI  *  •  •  #  O 

e  »\  *<w  *o  t  • 

—  r.  ircc  *c  •  ►■ 

•  mu  <v  •  4  •  •  •• 

►  •  *  •  •  O  tf'  ft.  •  • 

•M  M  Q  4*0*  *•  • 

—  C  »  •  i  N4 

•  •  V  •  •«.•_>  •  •  4 

—  in  *  »  »ipni  •  •  — 

o  •  n  c  trc  ••  c  c  • 

*•  r-ir  >  c  •  4  •  ►  i 

—  m  *,  r  <  «  •  o  ir*  m  —  •  v 

M  ••  .«  N  I  •  -  •  S. 

X  «  4  4  N  •  •  •  •  N  «. 

•  n  4  9>  •  >  N  •  W 
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